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ABSTRACT 


Analog-To-Digital  Converters  (ADCs)  are  integral  building  blocks  of  most  sensor  and 
communication  systems  today.  As  the  need  for  ADCs  with  faster  conversion  speeds  and 
lower  power  dissipation  increases,  there  is  a  growing  motivation  to  reduce  the  number  of 
power-consuming  components  by  employing  folding  circuits  to  fold  the  input  analog 
signal  symmetrically  prior  to  quantization  by  high-speed  comparators.  These  properties 
of  low-power  consumption,  compactness,  high-resolution  and  fast  conversion  speeds 
make  folding  ADCs  an  attractive  concept  to  be  used  for  defense  applications,  such  as 
unmanned  systems,  direction-finding  antenna  architectures  and  system-on-a-chip 
applications. 

In  this  thesis,  a  prototype  of  an  optical  folding  ADC  was  implemented  using  the 
Robust  Symmetrical  Number  System  (RSNS).  The  architecture  employs  a  three-modulus 
(Moduli  7,  8,  9)  scheme  to  preprocess  the  antenna  signal. 

This  thesis  focuses  on  the  simulation  and  hardware  implementation  of  this  ADC 
architecture,  including  the  bank  of  comparators  and  the  RSNS-to-Binary  Conversion 
within  a  Field  Programmable  Gate  Array  (FPGA),  to  achieve  an  eight-bit  dynamic  range 
of  133.  This  is  then  integrated  with  the  front-end  photonics  implementation  (designed 
under  a  separate  thesis). 

Low  frequency  analyses  of  the  results  using  a  1-kHz  input  signal  indicate  a  5.39 
Effective  Number  of  Bits  (ENOB),  a  Signal-to-Noise  Ratio  plus  Distortion  (SINAD)  of 
34.21  dB,  and  a  Total  Hannonic  Distortion  (THD)  of  -61.68  dB. 
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EXECUTIVE  SUMMARY 


Analog-To-Digital  Converters  (ADCs)  are  integral  building  blocks  of  most  sensor  and 
communication  systems  today.  As  the  need  for  ADCs  with  faster  conversion  speeds  and 
lower  power  dissipation  increases,  there  is  a  growing  motivation  to  reduce  the  number  of 
power-consuming  components  by  employing  folding  circuits  to  fold  the  input  analog 
signal  symmetrically  prior  to  quantization  by  high-speed  comparators.  The  folding  of  the 
analog  signal  allows  comparators  to  be  repetitively  used,  resulting  in  a  smaller  die  area 
and  lower  power  consumption. 

These  properties  of  low-power  consumption,  compactness,  high-resolution  and 
fast  conversion  speeds  make  folding  ADCs  an  attractive  concept  to  be  used  for  defense 
applications  that  involve  power  and  size  constraints  as  key  factors  in  the  design  of 
battlefield  systems  and  sensors. 

In  this  thesis,  a  prototype  of  an  optical  folding  ADC  was  implemented  using  the 
Robust  Symmetrical  Number  System  (RSNS),  which  minimizes  the  number  of 
comparators  and  removes  the  interpolation  circuitry  completely.  The  architecture 
employs  a  three-modulus  (Moduli  7,  8,  9)  scheme  to  preprocess  the  antenna  signal  and  is 
shown  in  Figure  1. 


Figure  1.  Block  Diagram  of  a  Three-Channel  Folding  ADC  Architecture. 
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The  goal  of  this  thesis  is  to  conduct  hardware  and  software  implementation  of  the 
Digital  Decoding  Sub-System  (DDS)  module  of  the  folding  ADC  architecture  (for 
Moduli  7,  8,  9),  from  the  bank  of  comparators  to  the  RSNS-to-binary  conversion  within 
the  Field  Programmable  Gate  Array  (FPGA),  as  well  as  integration  with  the  front-end 
Photonics  Encoding  Sub-System  (PES)  module  of  this  ADC  design.  This  was 
accomplished  via  several  milestones  described  below. 

Firstly,  the  RSNS  Dynamic  Range  (DR)  Computation  algorithm  was  verified  to 
be  correct,  proving  that  an  eight-bit  DR  of  133  can  be  achieved  theoretically  for  a  three- 
channel  RSNS  ADC  with  Moduli  m \  =  7,  m2  =  8  and  m 3  =  9. 

Secondly,  the  RSNS-to-binary  conversion  algorithm  for  a  three-channel  RSNS 
ADC  with  Moduli  m\  =  7,  m2  =  8  and  /h3  =  9  was  developed  in  Lab  VIEW,  shown  in 
Figure  2. 


UmbmouIk  Coda  PoadcaBit  CwdiboMi  Bit  E  Ttnitaiu*  AcTtr  Oilpa& 

ft'twn  ctmjirat-xj)  Cx.Ttrjwfi  Reveral  -  Dijpity 


Figure  2.  Lab  VIEW  Schematics  of  RSNS-to-Binary  Converter. 
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This  conversion  was  done  to  convert  the  RSNS  output  into  a  more  convenient 
decimal  representation.  This  implementation  was  shown  to  achieve  the  DR  value  of  133 
with  no  ambiguities,  which  is  in  agreement  with  the  RSNS  DR  Computation  algorithm, 
and  a  one-bit  improvement  over  that  achieved  in  a  previous  design.  Design  of 
thermometer  code  generator  circuits  and  simulation  of  this  algorithm  were  carried  out  to 
verify  that  it  is  working  properly  before  connecting  to  actual  signals. 

Thirdly,  the  comparator  circuits  and  RSNS-to-binary  conversion  algorithm  were 
designed  and  implemented  on  the  FPGA  to  form  the  DDS  module,  shown  in  Figure  3. 


New  Focus 
Photo-detectors 
(Quantity:  3) 


NI-9215  Analog 

Input  Module 


Data  Rate:  100  kS/s 
Channels:  4 


A 


B 


C 


NI-9111  Chassis 

FPGA:  XILINX  Virtex-5  LX30  (550  MHz) 


FPGA  Output 


Figure  3.  DDS  Module  Setup 


The  comparator  speed  was  a  limiting  factor  to  the  ADC  system  speed  as  actual 
comparator  ICs  were  used  for  sampling  outside  the  FPGA  in  a  previous  design.  The  key 
improvement  made  to  the  DDS  module  setup  in  Figure  3  is  that  the  comparator  circuits 
now  reside  in  the  FPGA,  allowing  them  to  sample  at  a  rate  equal  to  the  FPGA  speed. 
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The  comparator  circuit  and  RSNS-to-binary  conversion  logics  were  also  ran  on 
the  FPGA  to  guarantee  a  higher  FPGA  execution  speed,  as  opposed  to  running  it  on  a 
National  Instruments  Real-Time  Controller  module  with  a  lower  processing  speed.  These 
two  factors  allow  the  ADC  to  achieve  an  overall  higher  sampling  frequency. 

Lastly,  the  DDS  module  was  integrated  with  the  front-end  PES  module  to  form  a 
folding  ADC  system  and  characterization  of  the  ADC  performance  was  carried  out. 
Analysis  of  the  results  attributed  the  dominant  noise  source  in  the  ADC  system  to 
quantization  noise,  with  the  ADC  remaining  resilient  to  errors  caused  by  other  additive 
noise  sources  and  comparator  sampling. 

This  electro-optic  RSNS  ADC  system  has  been  demonstrated  to  work  and 
produces  an  eight-bit  output  with  relatively  simple  hardware  and  software.  Due  to  the 
reduced  number  of  hardware  and  software  components,  the  energy  and  size  savings  make 
this  folding  ADC  design  appealing  for  defense  applications,  such  as  unmanned  systems, 
direction-finding  antenna  architectures  and  electronic  warfare  system-on-a-chip 
applications. 


LIST  OF  ACRONYMS  AND  ABBREVIATIONS 


ADC 

Analog  to  Digital  Converter 

ASIC 

Application  Specific  Integrated  Circuit 

CRIO 

Compact  Reconfigurable  Input/Output 

CRT 

Chinese  Remainder  Theorem 

DDS 

Digital  Decoding  Sub-System 

DMA 

Dynamic  Memory  Allocation 

DNL 

Differential  Non-Linearity 

DR 

Dynamic  Range 

ENOB 

Effective  Number  of  Bits 

FIFO 

First-In-First-Out 

FPGA 

Field  Programmable  Gate  Array 

IC 

Integrated  Circuit 

INL 

Integral  Non-Linearity 

I/O 

Input/Output 

LPS 

Least  Positive  Solution 

LSB 

Least  Significant  Bit 

NI 

National  Instruments 

PES 

Photonics  Encoding  Sub-System 

PRP 

Pair-Wise  Relatively  Prime 

RNS 

Residue  Number  System 

RSNS 

Robust  Symmetrical  Number  System 

RTC 

Real-Time  Controller 

XIX 


sc 


Sub-Channel 


SINAD  Signal-to-Noise  Ratio  plus  Distortion 

SNR  Signal-to-Noise  Ratio 

THD  Total  Harmonics  Distortion 

VI  Virtual  Instrument 

M  Longest  series  of  consecutive  non-ambiguous  RSNS  vectors 

Pf  Fundamental  period 

XRg  Vector  of  RNS  residues  at  position  g 

Xh  Vector  of  RSNS  residues  at  position  h 


Table  1.  List  of  Variables. 


Variable 

Description 

Subscript  i 

Subscript  k 

St 

Encoder  Output 

Bit  Order 

hi 

Adder  Output 

Bit  Order 

hk 

RNS  States 

Channel  Number 

Bit  Order 

mi 

Modulus 

Channel  Number 

>h 

Alignment  Logic  Output  /  LPS  Solution 

Bit  Order  / 

Gate  Number 

Pik 

Position  Bit 

Channel  Number 

Bit  Order 

sik 

Thermometer  Code  /  RSNS  States 

Channel  Number 

Bit  Order 

h 

Shift  Value 

Channel  Number 

xx 
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I.  INTRODUCTION 


A.  FOLDING-TYPE  ANALOG-TO-DIGITAL  CONVERTERS 

Analog-To-Digital  Converters  (ADCs)  are  integral  building  blocks  of  most  sensor 
and  communication  systems  today.  They  allow  analog  data  measured  in  the  real  world  to 
be  sampled  and  converted  into  quantized  levels  for  high-speed  digital  processing. 

As  the  need  for  ADCs  with  faster  conversion  speeds  and  lower  power  dissipation 
increases,  there  is  a  growing  motivation  to  reduce  the  number  of  power-consuming 
components  by  employing  folding  circuits  to  fold  the  input  analog  signal  symmetrically, 
prior  to  quantization  by  high-speed  comparators.  [1] 

The  folding  of  the  analog  signal  allows  comparators  to  be  repetitively  used, 
resulting  in  a  smaller  die  area  and  lower  power  consumption.  One  folding  technique 
employed  is  the  Robust  Symmetrical  Number  System  (RSNS),  which  minimizes  the 
number  of  comparators  and  removes  the  interpolation  circuitry  completely.  [1] 

These  properties  of  low-power  consumption,  compactness,  high-resolution  and 
fast  conversion  speeds  make  folding  ADCs  an  attractive  concept  to  be  used  for  defense 
applications,  such  as  unmanned  systems,  direction-finding  antenna  architectures  and 
electronic  warfare  system-on-a-chip  applications. 

B.  SUMMARY  OF  RECENT  RESEARCH 

It  was  demonstrated  in  [2]  that  a  three-channel  folding  ADC  of  Moduli  3,  4  and  5 
can  be  designed  using  an  efficient  pipelined  RSNS-to-binary  algorithm  to  produce  a  six- 
bit  Dynamic  Range  (DR)  of  43  while  utilizing  significantly  less  electronic  components 
than  other  equivalent  six-bit  ADC  architecture  designs.  The  implementation  of  this  ADC 
architecture  was  verified  to  have  the  same  DR  in  [3],  shown  in  Figure  1. 
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Figure  1.  Block  Diagram  of  a  Three-Channel  Folding  ADC  Architecture  (After  [3]). 


In  [4],  it  was  demonstrated  that  the  same  algorithm  can  be  extended  and  applied 
to  a  folding  ADC  of  higher  moduli  (Moduli  7,  8  and  9).  A  seven-bit  DR  of  126  was 
achieved  in  [4],  but  the  full  DR  of  133  could  not  be  achieved  due  to  ambiguities  in  the 
ADC  output.  The  next  step  is  to  solve  for  the  ambiguities  and  carry  out  an 
implementation  of  this  ADC  design. 

This  thesis  is  done  in  conjunction  with  another  Master’s  thesis  to  implement  a 
three-channel  electro-optical  folding  ADC  architecture  (for  Moduli  7,  8,  9)  [5].  It  focuses 
on  the  implementation  of  the  Digital  Decoding  Sub-System  (DDS),  which  involves  the 
hardware  and  software  implementation  of  the  bank  of  comparators  to  the  RSNS-to-binary 
conversion  within  a  Field  Programmable  Gate  Array  (FPGA),  to  achieve  an  eight-bit  DR 
of  133. 
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This  is  done  using  a  novel  solution  to  solve  for  the  ambiguities  in  the  ADC  output.  The 
DDS  module  is  then  integrated  with  the  front-end  Photonics  Encoding  Sub-System 
(PES),  which  is  detailed  in  [5]. 

C.  PRINCIPAL  CONTRIBUTIONS 

Given  the  works  summarized  in  the  previous  section,  the  principal  contributions 
of  the  research  in  this  thesis  are  three-fold. 

Firstly,  this  thesis  provides  a  derivation  of  the  theoretical  DR  and  its  start  and  end 
positions  for  a  folding  ADC  of  Moduli  7,  8  and  9.  This  is  based  on  the  RSNS  dynamic 
range  computation  algorithm  in  [6]  and  [7]  and  is  used  as  a  yardstick  comparison  to 
verify  with  actual  experimental  values  obtained. 

The  second  contribution  is  the  implementation  of  the  RSNS-to-binary  algorithm 
in  Lab  VIEW  to  achieve  an  eight-bit  DR  of  133,  which  is  a  one -bit  improvement  over  that 
achieved  in  [4].  Simulation  of  this  algorithm  is  carried  out  to  verify  that  it  is  working 
properly  before  connecting  to  actual  signals.  This  is  done  by  designing  test  circuits  to 
produce  the  three-channel  thermometer  codes,  which  are  then  supplied  to  the  algorithm  to 
obtain  the  simulated  DR  and  its  position. 

Thirdly,  this  thesis  documents  the  DDS  implementation  of  the  comparator  circuit 
design  and  RSNS-to-binary  conversion  algorithm  onboard  a  National  Instruments  (NI) 
FPGA  so  as  to  achieve  a  higher  sampling  frequency.  The  comparator  circuit  is  redesigned 
in  Lab  VIEW  to  allow  it  to  achieve  a  higher  speed  to  match  that  of  the  FPGA,  as  opposed 
to  using  actual  comparator  Integrated  Circuits  (ICs)  as  in  [3],  The  DDS  module  is  then 
integrated  with  the  front-end  PES  module  in  [5]  to  form  an  overall  folding  ADC 
architecture  and  tested.  Lastly,  an  analysis  and  characterization  of  the  ADC  performance 
is  carried  out. 

D.  THESIS  OUTLINE 

Having  covered  the  recent  research  efforts  for  a  folding  ADC  architecture,  we 
will  provide  the  reader  with  a  summary  of  the  RSNS  structure  and  principles  behind  the 
RSNS-to-binary  conversion  algorithm  in  Chapter  II.  A  good  understanding  of  these  two 

3 


areas  is  required  in  order  to  translate  these  into  actual  hardware  and  software 
implementation.  A  theoretical  derivation  of  the  DR  and  its  start  and  end  positions  for  a 
folding  ADC  of  Moduli  7,  8  and  9,  based  on  the  RSNS  dynamic  range  computation 
algorithm  in  [6]  and  [7]  is  also  provided  in  Chapter  II. 

All  the  logic  blocks  and  their  equivalent  Boolean  equations  required  to  form  the 
RSNS-to-binary  converter  for  an  ADC  of  Moduli  7,  8,  9  are  delineated  in  Chapter  III. 

The  process  to  convert  the  Boolean  equations  in  Chapter  III  into  a  suitable  form 
for  implementation  in  Lab  VIEW  is  outlined  in  Chapter  III.  It  describes  how  each  logic 
block  is  designed  in  Lab  VIEW  and  the  integration  of  all  logic  blocks  to  form  the  RSNS- 
to-binary  converter. 

The  simulation  and  testing  of  the  RSNS-to-binary  converter  to  verify  that  the 
conversion  logic  is  functioning  properly  is  delineated  in  Chapter  V.  The  first  section  of 
this  chapter  describes  the  creation  of  thermometer  codes  for  each  channel  in  Lab  VIEW  to 
simulate  as  inputs  to  the  conversion  algorithm.  The  second  section  highlights  the  key 
simulation  results  obtained  for  various  models. 

The  process  of  implementing  the  comparator  circuit  design  and  RSNS-to-binary 
conversion  algorithm  on  a  NI  FPGA  module  is  detailed  in  Chapter  VI.  It  also  highlights 
the  integration  of  the  DDS  module  with  the  front-end  PES  module  in  [5]  to  form  an 
overall  folding  ADC  architecture  are  also  highlighted. 

An  analysis  of  the  results  to  characterize  the  performance  of  this  ADC  system  is 
provided  in  Chapter  VII. 

The  key  conclusions  obtained  from  this  project  and  recommendations  for  future 
research  are  given  in  Chapter  VIII. 
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II.  ROBUST  SYMMETRICAL  NUMBER  SYSTEM 


The  basic  theory  and  structure  of  the  RSNS,  and  how  its  properties  are  used  in  the 
implementation  of  the  RSNS-to-binary  converter,  is  explained  in  this  chapter. 
Equation  Section  (Next) 


A.  RSNS  STRUCTURE 

A  single-channel  RSNS  is  based  on  the  following  staircase  sequence  [1]: 

=  [0,1,2 -1  ,mi,mi  —  1, ..., 2, 1]  Equation  Section  (Next)  (1) 

This  sequence  starts  from  zero  and  increases  to  a  peak  value  m„  which  is  the  channel 
modulus.  It  then  decreases  back  to  zero  and  repeats  itself,  fonning  a  periodic  sequence 
with  a  period  of  2m,. 

A  three-channel  RSNS  vector  is  denoted  as: 


X,. 


(2) 


where  the  RSNS  residues  are  in  the  range 

-si  ={0, 1,  •••,  mj, 

s2  ={0, 1,  •••,m2},  (3) 

S3  ={0,1,  •  •  • ,  m3 } . 

For  an  A-channcl  RSNS,  each  number  in  the  sequence  is  repeated  N  times, 
extending  the  sequence  period  to  INnu.  Each  RSNS  channel  modulus  is  required  to  be 
pair-wise  relatively  prime  (PRP)  to  all  other  channel  moduli  [1]. 

The  three-channel  RSNS  case  (N=  3),  with  channel  moduli  mi  =  7,  m2  =  8  and  m3 
=  9,  which  meets  the  PRP  condition,  is  focused  on  in  this  thesis.  Each  channel  period  is 
of  length  2Nm,.  The  first  49  RSNS  vectors  of  the  three-channel  RSNS  structure  for  the 
moduli  mi  —  [7  8  9]  are  shown  in  Figure  2  as  an  illustration. 
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Figure  2. 


Three-Channel  RSNS  for  Moduli  mj  =  [7  8  9]. 


A  plot  of  the  folded  RSNS  waveforms  for  the  three  channels  in  Figure  2,  as  well 
as  the  relative  left-shift  for  Channels  2  and  3,  is  shown  in  Figure  3.  Notice  that  these 
relative  shifts  are  required  for  the  system  to  exhibit  Gray-code  properties,  where  the 
residues  within  consecutive  RSNS  vectors  change  one  at  a  time  at  the  next  code  position. 
This  property  makes  it  attractive  for  error  control  [1]. 


Channel  1  (m1  =  7) 


Figure  3.  Folded  RSNS  Waveforms  for  Moduli mi  =  [7  8  9]. 
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The  fundamental  period  of  a  three-channel  RSNS  is  [1]: 

P,  =  2N\{  mi  =  2Nm{m2m3  (4) 

This  is  calculated  tobePf  =  2(3)(7)(8)(9)  =  3024  for  a  RSNS  with mi  -  [7  8  9] . 

1.  Dynamic  Range 

The  maximum  DR  of  the  RSNS  is  defined  as  the  longest  series  of  consecutive 
unambiguous  RSNS  vectors  within  the  fundamental  period  of  the  system  [1],  The  DR  for 
the  same  RSNS  system  is  given  by  [1]: 

if  =  -m?  +  —  m,  +7.  (5) 

2  1  2  1 

'J  1  ^ 

This  is  calculated  as  if  =  —  (72)  +  —  (7)  +  7  =  133 . 

It  was  shown  that  the  size  of  the  DR  in  a  RSNS  is  the  same  regardless  of  the  shift 
sequence  and  whether  the  shift  is  to  the  left  or  right.  However,  the  location  of  the  DR  is 
affected  by  the  choice  of  the  channels  that  receive  the  shifts  [1]. 

Since  the  aim  is  to  achieve  an  eight-bit  DR  for  this  design  regardless  of  its 
location,  a  left-shift  tt  =  [0  1  2]  system  is  chosen  for  ease  of  implementation  as  the 
RSNS-to-binary  conversion  algorithm  is  derived  based  on  this  shift  system  in  [2],  with 
Channel  2  and  3  shifted  one  and  two  positions  to  the  left,  respectively,  both  relative  to 
Channel  1. 

2.  Ambiguity  Types 

The  conventional  approach  towards  finding  the  DR  and  its  position  is  to  search 
the  entire  fundamental  period  for  a  sequence  of  non-repeating  or  unambiguous  vectors  of 
unknown  length,  which  is  computationally-intensive.  A  more  efficient  approach  is  to 
compute  the  finite  locations  of  the  ambiguous  RSNS  vectors  for  each  channel  and  solve 
for  all  vector  ambiguity  locations  to  obtain  the  DR  and  its  position  [2],  which  will  be 
elaborated  in  Section  B  of  this  chapter. 
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Three  types  of  ambiguities  exist  in  each  RSNS  channel  [2].  The  positions  of  the 
three  ambiguity  types  for  the  Channel  1  (mi  =  7)  case  are  illustrated  in  Figure  4.  Type  0 
ambiguities  occur  for  every  repeating  channel  period,  while  Type  1  ambiguities  occur  on 
the  rise  and  fall  of  a  channel  period.  In  addition,  Type  2  ambiguities  occur  every  time 
each  residue  value  is  repeated  three  times  within  the  period. 


Type  0  Ambiguity 


Figure  4.  Single  Channel  Ambiguity  Types  ( ni\  =  7). 


3.  Sub-Channel  Analysis 

It  was  shown  in  [2]  that  decimating  each  of  the  three  channels  into  their  sub¬ 
channels  aids  in  the  solving  of  channel  ambiguities.  This  is  illustrated  for 
the  mt  -  [7  8  9]  case  in  Figure  5. 


Grouping  each  of  the  sub-channels  together,  and  re-indexing  the  position  index  h  to 
a  new  index  g,  we  obtain  the  results  in  Figure  6.  The  relationship  between  the  old 
position  index  h  and  the  new  index  g  is  given  by  [2]: 


g  =  - 


h 


g 

g  = 


3 

(h~  2) 


(Sub-Channel  0) 
(Sub-Channel  1) 
(Sub-Channel  2) 


(6) 
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Figure  5.  Decimation  of  Channels  into  Sub-Channels  for  Moduli  mi  =  [7  8  9]. 
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A  plot  of  the  RSNS  sub-channel  waveforms  in  Figure  6  is  shown  in  Figure  7. 
Each  channel  has  N  =  3  sub-channels.  From  Figure  7,  it  can  be  seen  that  Type  2 
ambiguities  are  eliminated,  leaving  only  Type  0  and  1  ambiguities. 

Sub-Channel  0 
Sub-Channel  1  _  _  . 


Channel  1  (/;/,  =  7)  Sub-Channel  2 


Figure  7.  Plot  of  Sub-Channels  0,  1  and  2. 


As  the  waveforms  are  folded  and  symmetric,  each  residue  value  occurs  twice  in  a 
single  folding  period,  except  for  the  minimum  and  the  maximum  value,  which  occur  only 
once.  Using  this  fact,  we  obtain  the  congruence  equations  describing  the  position  of  a 
RSNS  residue  vector  for  Sub-Channel  0  as  [2]: 

g  =  .v,  (mod  2m1 )  or  g  =  2m1  -  sx  (mod  2 ml ), 

g  =  s2  (mod  2 m2)  or  g  =  2m2  - s2  (mod  2 m2),  (7) 

g  =  s3  (mod  2m, )  or  g  =  2m3-s3  (mod  2 in, ). 
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The  congruence  equations  for  Sub-Channel  1  are: 

g  =sx  (mod  2 ml )  or  g  =  2 ml  -  sl  (mod  2m, ), 

g  =  s2  (mod  2 m2 )  or  g  =  2 m2  -  s2  (mod  2 m2 ),  (8) 

g  =  s3  -  I  (mod  2m, )  or  g  =  2m,  -s3  -1  (mod  2m, ). 

The  congruence  equations  for  Sub-Channel  2  are: 

g  =  .s',  (mod  2m, )  or  g  =  2m,  -  .v,  (mod  2m, ), 

g  =  s2  -  I  (mod  2 m2)  or  g  =  2m2  -j2  -1  (mod  2m2),  (9) 

gs^-1  (mod  2 m3)  or  g  =  2 m3  - .s'3  - 1  (mod  2m3). 

Equations  (7)-(9)  show  that  there  are  three  equations  for  each  sub-channel,  with 

■J 

two  choices  for  each  equation.  Thus,  each  RSNS  residue  vector  can  produce  up  to  2  or 
eight  unique  systems  of  equations.  This  means  that  any  RSNS  vector  can  have  up  to  eight 
ambiguities  within  the  fundamental  period. 

The  equations  also  reveal  that  the  start  position  of  the  DR,  h,  can  be  found  by 
detennining  the  sub-channel  of  a  particular  RSNS  vector,  solving  the  equivalent  set  of 
equations  to  obtain  g,  and  converting  g  to  find  h,  using  (6).  To  determine  which  sub¬ 
channel  a  RSNS  vector  is  from,  the  even-odd  structure  of  each  sub-channel  has  to  be 
investigated. 

4.  Even-Odd  Analysis 

The  even-odd  structures  ( e  =  even,  o  =  odd)  of  each  sub-channel,  and  the  overall 
even-odd  structure  of  a  three-channel  RSNS  withm,  =  [7  8  9] ,  are  illustrated  in  Figures 
8  and  9. 

From  Figure  8,  it  can  be  seen  that  each  sub-channel  produces  RSNS  vectors  with 
two  unique  even-odd  structures.  Thus,  the  Sub-Channel  0  equations  in  (7)  can  be  applied 

if  the  RSNS  vector  is  found  to  be  of  form  [ e  e  e]?  or  [o  o  o]r .  Similarly,  the  Sub-Channel 
1  equations  in  (8)  are  used  if  the  RSNS  vector  is  of  form  [e  e  o]r  or  [o  o  e]r  ,  and  the  Sub- 
Channel  2  equations  in  (9)  are  applied  if  the  RSNS  vector  is  of  form  [e  o  o]r  or  [o  e  e ]r  . 
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Figure  8.  Even-Odd  Structure  of  RSNS  Vectors  for  Sub-Channels  0,  1  and  2. 


This  even-odd  structure  repeats  in  blocks  of  six,  after  cycling  through  the  set  of 
sub-channels  twice,  as  shown  in  Figure  9.  This  means  that  the  minimum  distance 
between  ambiguous  parity  vectors  is  always  a  multiple  of  six  for  a  three-channel  RSNS, 
and  a  multiple  of  2 N  for  the  V-channcl  case. 
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Figure  9.  Overall  Even-Odd  Structure  of  RSNS  Vectors. 


It  was  shown  in  [2]  that  for  the  Sub-Channel  1  and  2  cases,  the  equations  in  (6) 
and  (7)  actually  convert  the  even-odd  structure  of  the  RSNS  vectors  to  the  form  [ e  e  e]  ?  or 

[o  o  o]r .  Thus,  the  Sub-Channel  0  case  can  be  exploited  as  a  base  case  to  develop  an 
efficient  RSNS-binary  conversion,  explained  in  Section  C  of  this  chapter. 
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B.  RSNS  DYNAMIC  RANGE  SEARCH  ALGORITHM 

An  efficient  RSNS  dynamic  range  search  algorithm  was  developed  in  [6]  and  [7], 
which  makes  use  of  the  cyclical  and  symmetrical  properties  of  the  RSNS  structure  to  cut 
down  on  the  solution  space  when  searching  for  the  DR.  This  is  faster  than  searching 
through  the  whole  solution  space,  which  is  computationally  intensive.  A  theoretical 
derivation  of  the  DR,  including  its  start  and  end  positions,  for  a  folding  ADC  of  Moduli 
7,  8  and  9,  based  on  this  algorithm,  is  provided  in  this  section. 

1.  Dynamic  Range  Upper  Bound 

The  upper  limit  of  M  is  defined  as  [6]: 

M  =  N  min  TT m  +  2TT m •  -1  (10) 

v,  7_ 

The  DR  upper  bound  for  an  ADC  of  Moduli  7,  8  and  9  is  calculated  to  be  22 1 : 

"m]  =  3[l  +  2(7-8-9)]-l  =  3026 
M]  =3[7  +  2(8-  9)]  - 1  =  452 
’m]  =3[8  +  2(7-9)]  — 1  =  401 
’m]4=3[9  +  2(7-8)]-1  =  362 

’m]  =3[(7-8)  +  2(9)]  — 1  =  221  (11) 

~M~\  =  3[(7  -  9)  +  2(8)]  - 1  =  236 
M]  =3[(8-9)  +  2(7)]  — 1  =  257 
m"|  =  3[(7- 8-9)  +  2(1)]  —  1  =  1517 
.'.  M~|  =  min|~M~|  =  221 

2.  RSNS  Vector  Ambiguity  Locations 

The  solutions  to  the  RSNS  vector  ambiguity  locations  for  a  three-channel  case  are 
shown  in  Figure  10.  Note  that  all  of  the  ambiguities  smaller  than  the  fundamental  period 

are  symmetric  around  a  Center  of  Ambiguity  (COA),  as  intuitively  shown  in  Figure  4  [6]. 
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The  three  identifier  digits  in  the  Case  Label  refer  to  the  number  of  Type  1 
ambiguities,  combination  number  and  sub-channel  number  respectively.  They  are  used  to 
provide  a  logical  case  numbering  system  for  implementation  in  a  computer  algorithm  [2]. 


Case 

Label 

Ambiguities  occur  at  h  and  h+k ,  where 
h  is 

k  is  a 
multiple 
of 

COR 

010 

Any  position  in  the  fundamental  period 

6m  -1  m2m3 

N/A 

110 

h  =  a  (3m  -i)  -  k  12 

6  m  2m3 

a  (3m  -1) 

111 

h  =  a  (3m  1  -  k  12 

a  (3m  1)  +  1 

112 

h  =  a  (3m  +  2  -  k  12 

a  (3m  1)  +  2 

120 

h  =  a  (3m  2)  -  k  12 

6m  im3 

a  (3m  2) 

121 

h  =  a  (3m  2)  +  1  -  k  12 

a  (3m  2)  +  1 

122 

h  =  a  (3m  2)  -  1  -  k  12 

a  (3m  2)  -  1 

130 

h  =  a  (3m  3)  -  k  12 

6m  ,m2 

a  (3m  3) 

131 

h  =  a(3m3)  -  2  -  k  12 

a  (3m  3)  -  2 

132 

h  =  a(3m3)  -  1  -  k  12 

a  (3m  3)  -  1 

210 

h  -  a  (3m  -im2)  -  k  12 

6m  3 

a  (3m  -im2) 

211 

h  =  a (3m  ,m2)+  /rs1  -  k  12 

a(3m  im2)  +  /?s1 

212 

h  =  a(3m  im2)  +  hs2  -  k  12 

a  (3m  t  m  2)  +  hS2 

220 

h  -  a(3m-|m3)  -  /c/2 

6m  2 

a(3m  -\m3) 

221 

h  -  a (3m  im3)  +  /7s1  -  /c/2 

a  (3m  -im3)  +  /is1 

222 

h  -  a(3m  im3)  +  hs2  -  k  12 

a(3m-im3)  +  /7s2 

230 

h  =  a(3m2m3)  -  /c/2 

6m  t 

a  (3m  2m  3) 

231 

h  -  a(3m2m3)  +  /rs1  -  /c/2 

a(3m2m3)  +  /7s1 

232 

h  -  a(3m2m3)  +  hs2-  k/2 

a(3m2m3)  +  /7s2 

310 

h  -  a  (3m-i m  2m  3)  -  k  12 

6 

a(3m-im2m3) 

311 

h  =  a (3m  im2m3)  +  /7s1  -  k/2 

a  (3m  -1  m2m3)  +  /?s1 

312 

h  =  a  (3m  t  m2m3)  +  hs2  -  k  /2 

a  (3m  1m2m3)  +  /?s2 

Figure  10.  RSNS  Vector  Ambiguity  Locations  (Three-Channel  Case). 


From  Figure  10,  it  can  be  seen  that  the  ambiguities  in  an  RSNS  structure  occur  at 
position  h  and  h  +  k,  where  k  is  a  multiple  of  the  moduli  combinations  for  each  case. 
There  is  an  inverse  relation  between  the  spacing  between  the  ambiguous  vectors  (k)  and 
the  spacing  between  the  COA. 
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Recall  that  a  left-shift  shifti  =  [0  1  2]  system  was  implemented  to  exhibit  gray- 
code  properties.  Thus,  the  base  cases  involving  Sub-Channel  0  (Case  XXO)  has  a  COA 
shift  of  hs0=  0.  The  cases  involving  Sub-Channel  1  (Case  XXI)  and  Sub-Channel  2 

(Case  XX2)  will  have  a  COA  shifts  of  Ki  and  Ki,  respectively.  These  are  the  least 
positive  solutions  to  the  following  two  sets  of  congruence  equations  [7]: 


~~~  s  0  (mod  7) 

~~~  =  0  (mod  8) 
h  +2 

— — ^ —  =  0  (mod  9) 


=  0  (mod  7) 

h  +1 

— ^ —  =  0  (mod  8) 
h  +1 

— ^ —  =  0  (mod  9) 


(12) 


The  COA  shifts  are  computed  to  be  hsl  =  1  and  hs2  =  23  using  the  generalized 
Chinese  Remainder  Theorem  (CRT)  procedure  described  in  the  Appendix. 

The  solutions  to  the  RSNS  vector  ambiguity  locations  for  a  RSNS  with  Moduli  7, 
8  and  9,  are  derived  from  Figure  10  and  summarized  in  Figure  11.  The  rows  highlighted 

in  grey  have  ambiguity  pairs  with  a  length  greater  than  If  =221  and  can  be  ignored  in 
the  DR  computation. 
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Case 

Label 

Ambiguities  occur  at  h  and  h+k ,  where 
h  is 

Ac  is  a 
multiple 
of 

COA 

010 

Any  position  in  the  fundamental  period 

3024 

N/A 

110 

h  =  a(21)  -  216 

432 

a  (21) 

111 

h  =  a( 21)  -  215 

a  (21)  +  1 

112 

h  =  a  (21)  -  214 

a  (21)  +  2 

120 

h  =  a  (24)  -  189 

378 

a  (24) 

121 

h  =  a  (24)  -  188 

a  (24)  +  1 

122 

h  =  a  (24)  -  190 

a  (24)  -  1 

130 

h  =  a  (27)  -  168 

336 

a  (27) 

131 

h  =  a  (27)  -  170 

a  (27) -2 

132 

h  =  a  (27)  -  169 

a  (27)  -  1 

210 

h  =  a  (168)  -  27 

54 

a  (168) 

211 

h  =  a (168)  +  /?s1  -  27 

a(168)  +  bs1 

212 

h  =  a  (168)  +  h  s2  -  27 

a  (168)  +  hs2 

220 

h  =  a  (189)  -  24 

48 

a  (189) 

221 

h  =  a (189)  +  /?s1  -  24 

a(189)  +  bs1 

222 

h  =  a (189)  +  hs2  -  24 

a  (189)  +  hs2 

230 

h  =  a  (21 6)  -  21 

42 

a  (21 6) 

231 

h  =  a (21 6)  +  hs1  -  21 

a(216)  +  bs1 

232 

h  =  a  (21 6)  +  h  s2  -  21 

a  (21 6)  +  bs2 

310 

h  =  a(  1512)  -  3 

6 

a(1512) 

311 

h  =  a(  1512)  +  /?s1  -  3 

a(1512)  +  /?s1 

312 

h  =  a(1512)  +  hs2  -  3 

a(1512)  +  hs2 

Figure  11.  RSNS  Vector  Ambiguity  Locations  (for  Moduli  7,  8  and  9). 


3.  Minimal  Ambiguity  Pair  Locations 

Due  to  the  symmetry  of  the  RSNS  vector  ambiguity  locations  about  Pf  /  2  ,  only 
ambiguity  pairs  from  h  =  -N  to  {Pf  /  2 )  +  N  need  to  be  considered  when  computing  M 
[7].  All  the  minimal  ambiguity  pair  (h{,h2)  locations,  derived  by  substituting  integer 
values  of  a  ,  hsl  =  1  and  hs2  =  23  into  Figure  1 1  for  each  of  the  cases,  are  shown  in  Figure 
12,  for  h  =  - 3  to  1515. 
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Figure  12. 


Minimal  Pair  Locations. 


4.  Consecutive  Minimal  Pair  Locations 

In  Figure  12,  all  minimal  pairs  with  a  starting  position  hx  earlier  than  the  starting 
position  of  the  previous  pair  are  removed.  The  remaining  minimal  pairs  are  sorted  such 
that  h2  is  monotonically  increasing  and  are  defined  as  consecutive  minimal  pairs  [7]. 
These  consecutive  minimal  pairs  are  derived  from  Figure  12  and  shown  in  Figure  13. 
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All 

COA 

h  2 

Case 

Distance 

Between 

Consecutive 

Minimal  Pairs 

-3 

0 

3 

310 

N.A. 

4 

25 

46 

231 

48 

82 

106 

130 

221 

125 

83 

107 

131 

222 

48 

141 

168 

195 

210 

111 

142 

169 

196 

211 

54 

165 

189 

213 

220 

70 

194 

215 

236 

232 

70 

195 

216 

237 

230 

42 

220 

241 

262 

231 

66 

271 

295 

319 

221 

98 

272 

296 

320 

222 

48 

309 

336 
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210 
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310 
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364 

211 

54 

332 

359 
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212 

75 

354 

378 
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220 
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410 
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432 
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42 
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457 

478 

231 
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508 
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669 

230 

42 

670 

673 

676 

311 

48 

h  i 

COA 

h2 

Case 

Distance 

Between 
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732 
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780 
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860 
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866 
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232 
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42 
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231 

66 
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1371 

210 

70 

1318 

1345 

1372 

211 

54 

1340 

1367 

1394 

212 

75 

1405 

1429 

1453 

221 

112 

1406 

1430 

1454 

222 

48 

1509 

1512 

1515 

310 

108 

Figure  13.  Consecutive  Minimal  Pair  Locations. 


The  DR  %f  is  the  largest  distance  between  endpoints  of  two  consecutive  minimal 
pairs,  computed  in  Figure  13.  The  result  is  a  fi/  =  133  starting  at  /;  =  /?,  +  !  =  733  and 
ending  at  /?  =  /?,- 1  =  865  ,  which  is  in  agreement  with  (5). 
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c. 


RSNS-TO-BINARY  CONVERSION  ALGORITHM 


With  a  good  understanding  of  the  underlying  RSNS  structure,  an  efficient  RSNS- 
to-binary  conversion  can  be  achieved  by  exploiting  the  relationship  between  the  RSNS 
and  the  Residue  Number  System  (RNS)  and  using  the  RNS  Least  Positive  Solution  (LPS) 
and  positional  alignment  techniques  to  solve  for  the  DR  position.  These  are  summarized 
in  this  section. 

1.  RSNS-RNS  Relationship 

It  was  shown  in  [2]  that  there  is  a  one-to-one  correspondence  between  the  RSNS 
and  the  RNS.  This  is  the  key  to  achieving  an  efficient  RSNS-to-binary  conversion  as  the 
RNS  has  no  ambiguities  in  a  fundamental  period,  unlike  the  RSNS. 

This  one-to-one  correspondence  between  the  RSNS  and  RNS  vectors,  using  the 
Sub-Channel  0  case  as  a  base  case  for  conversion,  is  demonstrated  in  Figure  14.  Each 
RSNS  residue  is  converted  to  a  unique  RNS  residue  such  that  there  is  no  ambiguity 
within  a  single  channel  period. 
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Figure  14.  One-to-One  Correspondence  between  RSNS  and  RNS  Vectors. 


Another  useful  RNS  property  is  that  every  XRg  vector  is  either  all  even  or  odd, 
unlike  the  RSNS. 
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It  is  possible  to  simplify  the  conversion  by  considering  only  the  even  XRg  vectors 
(grey-colored  columns  in  Figure  14)  and  dividing  their  values  and  the  index  g  by  two. 
This  is  illustrated  in  Figure  15. 
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Figure  15.  RSNS  and  RNS  Vectors  for  Even  Index  g. 


From  Figure  15,  it  can  be  seen  that  the  RSNS  vectors  are  transformed  into  RNS 
vectors  with  the  same  PRP  moduli.  The  index  position  g/2  can  now  be  solved  directly 
using  the  standardized  CRT  [10].  However,  there  are  still  up  to  2  or  eight  systems  of 
equations  to  solve  in  order  to  find  the  positions  of  a  single  RSNS  vector  within  the 
fundamental  period.  This  is  due  to  the  symmetry  of  each  RSNS  channel  period  and  three 
sub-channel  structure  within  each  channel  period. 

This  process  can  be  simplified  by  limiting  the  solution  range  to  within  the  DR. 
The  problem  is  then  reduced  to  finding  the  least  positive  solution  of  the  eight  systems  of 
equations,  explained  in  the  next  section. 

2.  Dynamic  Range  Compression 

A  closed-fonn  solution  for  the  start  and  end  DR  positions  for  moduli  of  form  (2'  - 
1,  2',  2'  +  1)  was  derived  in  [6].  The  radix  r  is  3  for  an  RSNS  with  moduli  mi  =  [7  8  9]. 
The  start  and  end  positions  can  be  calculated  using  [6]: 

Start  Position  of  M:  hx  =3(23'"'  -2r_1  -  2'  )  +  l  (13) 

and 
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(14) 


End  Position  of M:  h2  =  3(23'"'  +  22'"1)  +  1. 


The  DR  range  is  from  hx  =  73  3  and  h2  -  8  6  5  for  a  RSNS  with 
mi  =  [7  8  9].  The  corresponding  RSNS  vectors  are  V733  =[6  4  if  and  V865  =[6  0  l]r 
using  the  method  described  in  [5]. 

Recall  that  this  algorithm  is  applied  to  the  [e  e  e]r  vectors  in  the  Sub-Channel  0 
case  only.  Thus,  the  closest  [ e  e  e]T  vectors  to  Xm  and  XS65  must  be  chosen  for  this 
conversion  to  work.  They  are  X732  =[6  4  8]r  and  X864  =[6  0  0]r  . 


At  this  juncture,  it  must  be  noted  that  the  chosen  start  position  hx  —  732  falls 

outside  the  DR,  which  will  produce  ambiguous  results  during  implementation.  This  is  a 
special  case  where  the  DR  does  not  start  with  a  Sub-Channel  0  vector  and  end  with  a 
Sub-Channel  2  vector,  as  required  in  Figure  9. 


There  are  two  methods  to  overcome  this.  The  first  is  to  truncate  the  DR  to  a 
multiple  of  six,  so  that  the  DR  will  always  start  with  a  Sub-Channel  0  vector. 
Alternatively,  the  boundary  (start  and  end)  vectors  have  to  be  treated  as  special  cases  in 
the  Lab  VIEW  implementation.  The  details  of  both  methods  will  be  addressed  in  Chapters 
IV  and  V. 

Applying  (6),  we  calculate  the  corresponding  start  and  end  indices  for  the  RNS 
XRg  vectors  as: 


g\ 


\ 

3 


g  2 


h 

3 


864 


=  288. 


(15) 


Using  the  results  from  Figure  14,  we  consider  only  the  even  RNS  XRg  vectors. 
The  corresponding  start  and  end  indices  for  the  RNS  XRgal  2  vectors  are  calculated  as: 


_  244  _  122  g2  _  288  _ 
2  2  ’  2  2 


(16) 
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Using  the  results  from  (16),  we  obtain  the  total  number  of  RNS  vectors  required 

to  find  the  LPS  for  the  index  g  as  +  1  =  1 44 -122 +  1  =  23,  shown  in  Figure 

16.  The  RNS  vectors  in  Figure  16  are  derived  from  Figure  15  by  extending  the 
g  7  2  index  from  122  to  144.  The  representation  of  the  DR  is  now  compressed  from  a 
length  of  133  vectors  to  23  vectors. 
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Figure  16.  RNS  Vectors  Spanning  DR  (Before  and  After  Shifting  Start  Position). 


This  compression  is  possible  because  using  only  Sub-Channel  0  vectors 
compresses  the  DR  by  three,  and  exploiting  its  even-odd  structure  yields  an  additional 
compression  factor  of  two.  Thus,  the  total  compression  factor  is  2 N,  and  the  number  of 
RNS  vectors  required  to  find  the  LPS  for  the  index  g  is  [2]: 


M 

“133' 

2  N 

6 

23. 


(17) 


22 


Another  useful  RNS  property  is  that  subtracting  a  particular  RNS  vector  from  all  other 
RNS  vectors  in  the  fundamental  period  shifts  the  start  of  the  RNS  sequence  to  the 
position  of  the  same  vector  [2]. 


From  this  fact,  the  RNS  vector  — —  =  — —  =  [3  2  51r  is  subtracted  from  all 

2  2 


XR  XR 

RNS  vectors  between  - —  and  - —  to  obtain  the  shifted  RNS  sequence  in  the 

2  2 


bottom  row  of  Figure  16.  This  means  that  the  solution  for  the  index  h  is  now  shifted  to 
the  range  0</?'<132  rather  than  733</?<865.  This  is  desirable  as  the  former  can  be 
represented  in  an  8-bit  binary  number,  while  the  latter  requires  10  bits. 


3.  Alignment  of  RNS  Least  Positive  Solution 

A  method  of  finding  the  LPS  to  multiple  systems  of  equations  is  a  positional 
alignment  solution.  This  method  asserts  the  positions  of  the  RNS  residues  for  each 
XR 

channel  of  the  — vectors,  shown  in  Figure  16.  The  LPS  is  the  position  g'/2of  the 
XR 

first  — vector  in  which  all  three  asserted  residues  align  [2], 

After  finding  the  LPS  g '/  2  for  a  particular  RSNS  vector  Xh  =  [.s’,  s2  ,s\]7  ,  the  LPS 
is  converted  back  to  the  shifted  index  h '  by  reversing  the  sub-channel  and  even-odd 
compression  process  carried  out  before  the  LPS  alignment.  Using  Fsc  to  denote  the  sub¬ 
channel  compensation  factor,  Fodd  as  the  even-odd  compensation  factor,  and  (6),  the 
expression  to  obtain  h  '  is  [2]: 

h '  =  3(g '+  Fodd)  +  Fsc,  (18) 

where 
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[  0,  if  residue  sl  from  Xh  is  even 
[  1,  if  residue  sl  from  Xh  is  odd 


F< 


sc 


0,  if  Xh  is  from  Sub-Channel  0 

1,  if  Xh  is  from  Sub-Channel  1  k 

2,  if  Xh  is  from  Sub-Channel  2 


(19) 


Since  h '  is  the  index  h  after  shifting  the  DR  start  position  to  zero,  the  position  of  the 
RSNS  vector  within  the  DR,  or  index  h,  can  simply  be  obtained  by  adding  the  DR  start 
position  h\\ 

h  =  hl+h\  (20) 

The  key  principles  and  steps  of  the  RSNS-to-binary  conversion  algorithm 
described  in  this  chapter  are  encapsulated  in  Figure  17.  An  efficient  RSNS-to-binary 
conversion  can  be  achieved  by  exploiting  the  RSNS-RNS  relationship,  as  well  as  using 
the  RNS  LPS  and  positional  alignment  techniques  to  solve  for  the  DR  position. 


Figure  17.  RSNS-to-Binary  Conversion  Algorithm  (After  [2]). 


The  translation  of  the  principles  of  this  RSNS-to-binary  conversion  algorithm  into 
a  feasible  logic  block  diagram  is  explained  in  the  next  chapter. 
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III.  RSNS-TO-BINARY  CONVERSION 


Given  the  key  principles  of  the  RSNS-to-binary  conversion  algorithm  described  in 
the  last  chapter,  the  translation  of  this  algorithm  into  a  feasible  logic  block  diagram  for 
implementation  in  the  Lab  VIEW  programming  environment  is  explained  in  this  chapter. 

A.  LOGIC  BLOCK  DIAGRAM  OF  RSNS-TO-BINARY  CONVERTER 

The  overall  logic  block  diagram  of  this  converter  that  can  be  implemented  in 
Lab  VIEW  and  the  input/output  variables  for  each  block  are  shown  in  Figure  18.  The 
subsequent  sections  will  detail  why  each  logic  block  is  required  and  how  it  is  derived. 


Figure  18.  Logic  Block  Diagram  of  RSNS-to-Binary  Converter  (After  [4]). 


The  inputs  of  this  conversion  system  come  from  a  hank  of  24  comparators  at  the 
end  of  a  photonic  analog  folding  circuit  and  can  be  represented  as  a  three-channel  RSNS 
vector  as  shown  in  (2).  These  comparator  outputs  are  separated  into  three  channels  of 
seven,  eight  and  nine  comparators,  according  to  the  Moduli  mj  =  [7  8  9]. 


25 


The  input  residues  [.s’l,.v2,.v3]7are  encoded  in  a  thermometer  code  for  this 
implementation,  according  to  (3).  Each  of  the  bits  in  the  thermometer  code  is  labeled 
as  sik  ,  where  the  index  i  represents  the  channel  of  the  RSNS  residue  vector,  and  the  index 
k  is  the  bit  position  in  the  thermometer  code,  with  k  —  0  corresponding  to  the  position  of 
the  Least  Significant  Bit  (LSB).  The  RSNS  thermometer  code  bits  for  mi  =  [7  8  9]  are 
shown  in  Figure  19. 
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Figure  19.  RSNS  Thermometer  Codes  for  mi  =  [7  8  9]. 
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B. 


POSITION  BIT  CONVERSION 


The  next  step  is  to  convert  the  RSNS  thermometer  code  residues  into  RNS 
residues  using  the  one-to-one  correspondence  property  of  the  RSNS-RNS  relationship. 
Each  RNS  residue  can  be  represented  by  a  unique  position  bit. 

The  position  bits  are  denoted  as  pjk ,  where  the  index  i  is  the  channel  and  k  is  the 
RNS  residue  value  with  0  <k<  ( mi  -1] .  The  LPS  is  then  found  by  finding  the  position 
g '/  2  where  the  three  position  bits  from  each  channel  are  asserted  and  aligned. 

1.  RSNS  Thermometer  Code  to  RNS  Residue/Position  Bit  Conversion 

The  simplest  case  for  this  conversion  is  to  consider  the  residues  for  Channel  1. 
From  Equations  (7),  (8)  and  (9),  we  see  that  the  top  row  of  each  set  of  equations  is 
identical.  This  means  that  the  conversion  of  RSNS  residues  to  RNS  residues  for  Channel 
1  is  the  same  regardless  of  the  sub-channel  that  the  residue  comes  from.  The  conversion 
process  for  Channel  1  ( ml  =  7)  is  shown  in  Figure  20. 
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Figure  20.  Channel  1  RSNS-to-RNS  Conversion  for  all  Sub-Channels  (After  [2]). 
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Recall  that  each  RSNS  residue  value  5,  can  be  represented  by  two  unique  RNS 
residues  to  prevent  ambiguity  within  a  single  channel  period,  as  illustrated  in  Figure  14 
previously.  This  conversion  is  shown  as  Step  0  to  1  in  Figure  20. 

The  RNS  residues  are  then  rounded  down  to  the  nearest  even  RNS  residues  and 
divided  by  two  to  obtain  RNS  residues  with  the  same  PRP  moduli  as  the  RSNS,  as 
illustrated  in  Figure  15  previously.  This  conversion  is  shown  as  Steps  1  to  3  in  Figure  20. 
The  asserted  position  bits  corresponding  to  each  of  the  RSNS  and  RNS  residues  are 
shown  in  Figure  2 1 . 
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Channel  1  RSNS-RNS-Position  Bit  Correspondences  for  all  Sub- 
Channels  (After  [2]). 


The  conversion  process  for  Channel  2  ( m2  =  8)  and  Channel  3  (m3  =  9)  are 
similar  to  that  of  Channel  1,  except  that  they  are  left-shifted  by  one  and  two  positions, 
respectively.  Examination  of  (7),  (8)  and  (9)  reveals  that  this  impacts  the  conversion  of 
RSNS  residues  from  different  sub-channels. 

The  RSNS  to  position  bit  conversion  process  for  Sub-Channels  0  and  1  of 
Channel  2  are  illustrated  in  Figures  22  and  23.  The  same  process  for  Sub-Channel  2  is 
illustrated  in  Figures  24  and  25. 
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Note  that  Sub-Channel  2  is  shifted  by  one  relative  to  Sub-Channels  0  and  1,  from 
(7),  (8)  and  (9).  A  comparison  of  Figures  23  and  25  shows  that  the  position  bits  of  Sub- 
Channel  2  are  the  reversed  of  those  for  Sub-Channels  0  and  1 . 
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Figure  22. 

Channel  2  RSNS-to-RNS  Conversion  for  Sub-Channels  0  and  1  (After 

[2]). 
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Figure  23. 
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Figure  24.  Channel  2  RSNS-to-RNS  Conversion  for  Sub-Channel  2  (After  [2]). 
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Figure  25.  Channel  2  RSNS-RNS-Position  Bit  Correspondences  for  Sub-Channel  2 

(After  [2]). 
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The  RSNS  to  position  bit  conversion  process  for  Sub-Channel  0  of  Channel  3  are 
illustrated  in  Figures  26  and  27.  The  same  process  for  Sub-Channels  1  and  2  is  illustrated 
in  Figures  28  and  29.  Note  that  Sub-Channels  1  and  2  are  shifted  by  one  relative  to  Sub- 
Channel  0,  from  (7),  (8)  and  (9).  A  comparison  of  Figures  27  and  29  shows  that  the 
position  bits  of  Sub-Channels  1  and  2  are  the  reversed  of  those  for  Sub-Channel  0. 
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Figure  26. 

Channel  3  RSNS-to-RNS  Conversion  for  Sub-Channel  0  (After  [2]). 
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Figure  27.  Channel  3  RSNS-RNS-Position  Bit  Correspondences  for  Sub-Channel  0 

(After  [2]). 


31 


RSNS 

Residue 


RNS 

Residue 


Even  RNS  Residue 


RNS  Residue  for  PRP 
Moduli 


even(s3  - 1)  even(l 7  —  s3 ) 


S3 

s3  —  1 

17 -s3 

even(s3) 

even(17  —  s3) 

2 

2 

0 

17 

17 

16 

16 

8 

8 

1 

0 

16 

0 

16 

0 

8 

2 

1 

15 

0 

14 

0 

7 

3 

2 

14 

2 

14 

1 

7 

4 

3 

13 

2 

12 

1 

6 

5 

4 

12 

4 

12 

2 

6 

6 

5 

11 

4 

10 

2 

5 

7 

6 

10 

6 

10 

3 

5 

8 

7 

9 

6 

8 

3 

4 

9 

8 

8 

8 

8 

4 

4 

Figure  28.  Channel  3  RSNS-to-RNS  Conversion  for  Sub-Channels  1  and  2  (After 

[2]). 
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Figure  29.  Channel  3  RSNS-RNS-Position  Bit  Correspondences  for  Sub-Channels  1 

and  2  (After  [2]). 
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2.  Position  Bit  Equations 

Using  logic  tables  and  Karnaugh  mapping,  we  can  show  that  the  general 
equations  of  the  position  bits  for  even  and  odd  moduli  are  [2]: 

P,0=Sil’ 

Pil  ‘hi  Si3  ’ 

Pi2  ~  Si3  Si5  ’ 


P  j  |  mi  ^ 


'h'(  in.  -4  ) 


^ i(tni -2)  ’ 


for  even  moduli,  and 


P i(m,-3)  ~  Si4Si6 ’ 
2)  =  Si2Si4’ 
P i(m(— 1)  —  SiOSi2’ 

PiO  =Sn, 

Pn  =  SnSi3, 

Pil  =  Si3Si5’ 


Pi(m ,-3)  _  Si4Si6’ 

Pi[m,- 2)  = 

Pi(m,- 1)  =  SiOSi2> 


(21) 


(22) 


for  odd  moduli. 
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Equations  (21)  and  (22)  are  used  to  generate  the  position  bits  as  follows: 


for  ml  =  7  ,  and 


for  m2=  8 ,  and 


for  m3=  9 . 


Pi 0  ~  Sll> 

Pi  i  ~  i  *^13 5 

P 12  —  S\3S\5’ 

Pl3  ~  S15’ 

P\4  =  S14S16’ 
Pl5  ~  S12S14’ 

P 16  —  510512’ 

P  20  =  ^21’ 

Pl\  =  ^21^23’ 
Pl2  =  ^23  S25  ’ 

P 23  ~  S25  S21  ’ 
Pl4  ~  S26  ’ 
P25  ~  S24S26’ 
Pl6  =  S22  S24  ’ 

Pll  ~  S20  S22  ’ 

P30  =«31. 

P3 1  =  531  ^33’ 
P32  =  ^33  ^35  ’ 
P 33  =  S35  S31  ’ 

P34  =  S31  ’ 

P  33  =  ^36^38’ 
P36  =  ^34  ^36  ’ 
P 31  =  S32  ^34  ’ 
P3S  =  S30  S32  ’ 


(23) 


(24) 


(25) 
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c. 


EVEN  RESIDUE  AND  SUB-CHANNEL  FLAGS 


From  Figures  20-29,  it  can  be  seen  that  only  the  position  bits  representing  even 
RSNS  residues  for  all  three  channels  require  bit  reversal,  depending  on  the  sub-channel 
they  come  from.  This  means  that  a  logic  block  is  required  in  the  converter  diagram  to 
check  for  even  residues  and  their  sub-channels. 

1.  Even  Residue  Flags 

The  letter  V  is  used  to  represent  an  even  residue  flag,  with  the  subscript  denoting 
the  channel.  Each  flag  will  be  asserted  when  the  RSNS  residue  is  even.  It  was  shown  in 
[2]  that  the  general  equations  of  the  even  residue  flags  are: 


Even  Moduli:  e.  =  si0  +  snsi2  +  si3si4  +--  +  s.t 


(m,— l) 


(26) 


and 


Odd  Moduli:  et  =  si0  +snsi2  +si3si4  h - h  s. 


(mi  —2)  ^  i(  mi  -l) 


(27) 


Equations  (26)  and  (27)  are  used  to  generate  the  even  residue  flags  as  follows: 


Channel  1 :  e,  =  slQ  +5n512  +  s13sl4  +  sl5sl6, 

Channel  2:  e2  =  s20  +  s2ls22  +  s23s24  +  s25s26  +  s21,  (28) 

Channel  3:  e3  =  s30  +  sns32  +  s33s34  +  s35s36  +  s31s3S. 

2.  Sub-Channel  Flags 

The  symbol  ‘SC’  is  used  to  represent  a  sub-channel  flag,  with  the  subscript 
denoting  the  channel.  Each  flag  is  used  to  determine  if  the  position  bits  of  each  RSNS 
residue  require  reversal,  depending  on  their  sub-channels.  With  the  sub-script  N  to  denote 
the  number  of  channels  in  the  system,  it  was  shown  in  [2]  that  the  general  equations  of 
the  sub-channel  flags  are: 
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SCNA  =el  ®e2, 

SC2  —6n_ 2  ®eN_i, 

=eN®eN- 1’ 

SC0  -eN®el, 

where  ®  denotes  an  XOR  operation. 

Equation  (29)  is  used  to  generate  the  sub-channel  flags  as  follows: 


SC0  =  e3  ©  ex , 

SCX  =  e,  ©  e, ,  (30) 

SC2  =e2  ©er 

D.  CONDITIONAL  BIT  REVERSAL 

Conditional  bit  reversal  of  the  position  bits  is  based  on  the  sub-channel  flags.  In 
order  to  maintain  proper  housekeeping  of  variable-naming,  all  position  bits  after 
inversion  have  an  additional  subscript  ‘o’. 

Based  on  Figures  20-29,  the  position  bits  of  Channel  1  are  never  reversed.  The 
position  bits  of  Channel  2  are  reversed  if  the  residue  is  from  Sub-Channel  2,  i.e.,  when 
SC2  is  asserted.  The  position  bits  of  Channel  3  are  reversed  if  the  residue  is  from  Sub- 

Channel  1  or  2,  i.e.,  when  SC0  is  asserted.  Note  that  SCX  is  not  required  for  this  shift 
sequence. 

This  reversal  procedure  can  be  accomplished  via  multiplexer  circuits  using  the 
sub-channel  flags  as  control  signals  and  will  be  shown  in  Chapter  IV. 
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E.  LEAST  POSITIVE  SOLUTION  ALIGNMENT 

After  determining  the  position  bits  for  each  channel,  the  position  bits  from  all 
three  channels  are  then  used  to  compute  the  LPS  of  the  positional  alignment.  One  useful 
property  is  to  make  use  of  the  one-to-one  correspondence  between  the  RNS  residues  and 
the  position  bits.  This  allows  the  RNS  residues  to  be  replaced  by  the  position  bits  in 
Figure  16.  This  is  shown  in  Figure  30. 

Recall  from  (15)  and  (16)  that  the  RNS  vector  — = — ^~  =  [3  2  5]r 

m  n  *1  J1 

corresponds  to  the  RSNS  vector  X732  =[6  4  8]  ,  which  is  of  form  [ e  e  ej  and  from  Sub- 
Channel  0.  From  Figure  21,  the  RSNS  Channel  1  residue  5,  =  6  corresponds  to  RNS 
residues  of  3  and  4,  which  are  asserted  and  highlighted  in  Figure  30.  Similarly,  from 
Figure  23,  the  RSNS  Channel  2  residue  s2  =  4  corresponds  to  RNS  residues  of  2  and  6, 

and  the  RSNS  Channel  3  residue  =  8  corresponds  to  RNS  residues  of  3  and  5  from 
Figure  27. 
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Figure  30.  LPS  Alignment  using  Position  Bits. 
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The  asserted  RNS  residues  for  the  three  channels  are  shown  in  Figure  30.  Ideally, 
only  one  LPS  should  be  active  within  the  DR  to  prevent  ambiguities.  However,  there  is 
more  than  one  LPS  ( LPSO  and  LPS8 )  asserted  for  this  case. 

There  is  an  ambiguity  in  this  case  because  the  DR  length  (133)  is  not  evenly 
divisible  by  six,  and  the  DR  starts  at  an  index  position  (733)  that  is  not  evenly  divisible 
by  six.  There  are  two  methods  to  overcome  this,  and  these  are  explained  below. 

1.  Full  Dynamic  Range  With  LPS  Priority 

If  the  DR  is  not  truncated  to  a  length  that  is  evenly  divisible  by  six,  it  is  possible 
that  two  LPS  equations  will  be  asserted  at  the  same  time  -  one  of  the  boundary  LPS 
equations  (i.e.,  LPSO  or  LPS22 )  and  LPSX.  If  the  full  DR  is  to  be  maintained,  the 
boundary  LPS  equation  should  always  be  ignored  in  favor  of  the  LPSX  solution. 

This  can  be  done  by  using  additional  logic  circuitry  to  check  for  such  cases  and 
giving  priority  to  the  LPSX  solution  over  the  boundary  LPS  equations.  The  details  of  this 
circuitry  are  elaborated  in  Chapter  V. 


2.  Truncated  Dynamic  Range 


The  DR  can  be  truncated  to  a  length  that  is  evenly  divisible  by  six,  i.e.,  126,  so 
that  it  starts  with  a  Sub-Channel  0  vector  and  ends  with  a  Sub-Channel  2  vector,  as 
required  in  Figure  9.  The  number  of  RNS  vectors  required  to  span  the  DR  is  now: 


M 

“126“ 

2  N 

6 

(31) 


For  the  first  case  where  the  full  DR  is  maintained,  the  number  of  RNS  vectors 
required  to  span  the  DR  is  23,  in  accordance  with  (17).  The  full  23  alignment  equations 
are  given  by 
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LPS 0:  n0  =  pl3ap22ap35a, 

LPS 1:  nl=pl4ap23ap36a, 

LPS 2:  n2  =  pl5ap24ap37a, 

LPS 3:  n3  =  pl6ap25ap3Sa, 

LPS 4:  n4  =  pWap26ap30a, 

LPS 5:  n5  =  pUap27ap3la, 

LPS 6:  n6  =  pnap20ap32a, 

LPS1:  n7  =  pX3ap2Up33a, 

LPS 8:  n8=pl4ap22ap34a, 

LPS 9:  n9  =  pl5ap23ap35a, 

LPSlO:n10  =  Pl6ap24aP36a, 

LPSU:  nu  =  pl0ap25ap37a, 

LPS\2:nu=  pUap26apna, 

LPS  13:  n13  =  pnap27ap30a, 

LPSU:  nl4  =  pl3ap20ap3la, 

LPS15:  nl5  =  pUap2lap32a, 

LPS16:  nl6  =  pl5ap22ap33a, 

LPS\l:nxl  =  pl6ap23ap34a,  (32) 

ZP518:  «18  =  PWaP24aP35a, 

LPSU:  nX9  =  pnap25ap36a, 

LPS 20:  n20=p12ap26ap37a, 

LPS 21:  n21=pl3ap27ap38a, 

LPS22:  n22=pl4ap20ap30a. 

These  operations  can  be  accomplished  using  3-input  NAND  gates.  Note  that  the  last  two 

LPS  equations  ( LPS21  and  LPS22 )  are  not  required  for  the  truncated  DR  case. 
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F.  ENCODER 

The  next  step  is  to  convert  from  the  LPS  index  g  /  2  to  a  binary  representation  of 
the  position  h  of  the  RSNS  vector  within  the  DR  bounds.  Since  only  one  of  the  LPSX 
equations  will  be  asserted  at  any  one  time,  an  encoder  can  be  used  to  convert  the  23 
possible  outputs  to  a  binary  representation  using  B  bits,  calculated  via 

B  =  \og2\L\  (33) 

Five  bits  are  required  for  both  the  full  and  truncated  DR  cases,  with  a  23-to-5  encoder  for 
the  full  DR  case  and  a  21-to-5  encoder  for  the  truncated  DR  case. 

G.  ADDER 

As  the  output  of  the  encoder  is  a  five-bit  representation  of  g  /  2 ,  a  left-shifted 
version  of  g  /  2  is  equivalent  to  a  multiplication  by  two  and  converts  it  to  a  six-bit  index 
g.  From  (18),  the  Fodd  compensation  factor  has  to  be  computed  next,  as  the  DR  was 

compressed  using  even  residues  only.  The  factor  Fodd  is  equal  to  ex  since  the 
complement  of  ex  is  asserted  when  the  Channel  1  residue  5,  is  odd,  as  shown  in  (28).  The 
LSB  of  g  is  guaranteed  to  be  zero,  as  g  is  a  left-shifted  version  of  g  /  2 .  Thus,  the  even 
residue  flag  ex  can  replace  the  LSB  of  g. 

From  Equation  (18),  the  computation  of  3(g  +  e|)  is  also  required.  It  is  easier  to 
implement  this  as  3X  =  2X  +  X  in  hardware  using  a  wired  shift  and  an  adder  rather  than 
a  multiplication  by  three  using  a  multiplexer.  To  achieve  this,  (g  +  ej)  is  shifted  one 

position  left  to  form  2(g  +  e1)  and  used  as  the  input  to  an  eight-bit  carry-look- ahead 
adder.  This  concept  is  illustrated  in  Figure  31,  where  the  notation  gh  refers  to  the  bth  bit 
of  the  binary  representation  of  index  g  [2], 
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2(g  +  e i  ) 

S 4  §3 
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Si 

So  0 

_Xi+  -» 
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S2 

Si  So  ei  + 

3X 

2[g  +  ei  )  h 

7  ^6  ^5 

K 

h'2 

-  o 

‘-sT 

"-sT 

Figure  3 1 .  Adder  Function  for  Implementing  Multiplication  by  Three  (After  [2]). 

Lastly,  Equation  (18)  requires  the  addition  of  the  Fsc  compensation  factor,  as  the 
DR  was  compressed  using  Sub-Channel  0  vectors  only.  The  Fsc  compensation  factor  is 

represented  by  the  signals  SC 2  and  SO  from  (30),  which  are  asserted  if  the  Channel  2 
residue  is  from  Sub-Channel  2,  and  the  Channel  3  residue  is  from  Sub-Channel  1  or  2, 
respectively. 

Fortunately,  the  left  shift  of  (g  +  e,)  provides  a  LSB  slot  guaranteed  to  be  zero, 
and  the  carry-in  to  the  adder  provides  another  LSB  slot.  This  allows  (g  +  e J),  2(g  +  e1), 

SCO  ,  and  SC2  to  be  summed  in  a  single  adder  [2],  shown  in  Figure  32. 

2X  2(g  +  e,  )  g4  g3  g2  gj  g0  e,  SC o  <-  SC2  (carry-in) 

-2L+  (s  +  ei)  +  _ S  4  gj  §2  Si  go  e\  + 

3  A  3  ( g  +  ex  j  h7  h6  h5  h4  h3  h2  hl  h0 

Figure  32.  Single  Adder  for  Converting  LPS  to  Binary  (After  [2]). 

The  output  of  this  adder  is  the  binary  representation  of  the  position  h  within  the 
DR  for  the  RSNS  vector  Xh  -  [s,  s2  s3  ]' .  The  Lab  VIEW  circuit  schematics  that 
implement  the  logic  equations  developed  in  this  chapter  are  provided  in  the  next  chapter. 
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IV.  LABVIEW  IMPLEMENTATION  OF  RSN S-TO-BIN ARY 

CONVERSION 


Using  the  logic  equations  developed  in  the  previous  chapter,  we  document  the 
process  of  implementing  the  RSNS-to-binary  converter  in  the  NI  LabVIEW 
programming  environment  in  this  chapter. 

A.  LABVIEW  SCHEMATICS  OF  RSNS-TO-BINARY  CONVERTER 

The  NI  LabVIEW  programming  environment  was  specifically  selected  as  its 
programming  structure  is  modular  in  nature.  This  allows  the  codes  for  each  logic  block  to 
be  stored  as  sub-routines  or  sub- Virtual  Instruments  (sub- Vis)  and  run  as  part  of  a  larger 
routine  or  Virtual  Instrument  (VI). 

This  modular  structure  is  beneficial  for  the  future  expansion  of  such  RSNS  ADC 
system  as  it  allows  additional  sub-routines  to  be  designed  and  added  to  existing  routines 
should  it  be  decided  to  scale  the  ADC  to  higher  moduli  configurations  in  future. 

In  addition,  if  the  ADC  requires  future  upgrading,  it  is  relatively  easy  to  replace 
current  FPGA  modules  with  other  higher-speed  and  higher-bandwidth  NI  modules  with 
minimal  disruption  to  existing  codes  due  to  its  plug-and-play  features. 

The  overall  schematics  of  this  converter  in  LabVIEW  are  shown  in  Figure  33 
based  on  the  major  logic  blocks  developed  in  Chapter  III.  The  process  of  constructing 
each  logic  block  as  a  sub-routine  and  implementing  it  as  part  of  the  overall  routine  is 
detailed  in  the  subsequent  sections. 
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Figure  33.  Overall  Lab  VIEW  Schematics  of  RSNS-to-Binary  Converter. 


B.  POSITION  BIT  CONVERSION  AND  EVEN  RESIDUE/SUB-CHANNEL 

FLAGS 

The  logic  blocks  for  position  bit  conversion  and  even  residue  flags  are  combined, 
as  both  of  them  can  be  calculated  using  the  thermometer  codes  from  the  comparators. 
The  Boolean  equations  used  are  modified  from  those  in  Chapter  III.  The  only  changes  are 
the  application  of  DeMorgan’s  Theorem  to  use  NAND,  NOR  and  inverter  gates 
predominantly. 

1.  Channel  1 

DeMorgan’s  Theorem  was  applied  to  equations  (23)  and  (28)  to  obtain: 
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A  o 
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=  511 
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Pn 

=  %515  = 

=  % 

+ 
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Pn 

=  S15’ 

Pn 

=  susie  z 

=  *14 

+ 

^16 5 

Pn 

=  ^12  ^14  = 

=  512 

+ 

^14’ 

Pn 

=  ^10^12  = 

=  ‘ho 

+ 

S12> 

e\  510  +  511  SU  +  ^13  ^14  +  ^15  ^16  ^10 


(34) 


(35) 


The  implementation  of  equations  (34)  and  (35)  in  Lab  VIEW  is  shown  in  Figure  33. 


slO  sll  s!2  s!3  s!4  s!5  s!6 


Figure  34.  Lab  VIEW  Schematics  of  Channel  1  Position  Bit  and  Even  Residue  Flag. 
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2. 


Channel  2 


DeMorgan’s  Theorem  was  applied  to  equations  (24)  and  (28)  to  obtain 


and 


P  20 

~  S2l’ 

Pi  1 

=  ^21^23 

=  Vl 

4"  . 

S23> 

P  22 

=  S  23S  25 

=  S23 

+ 

S25’ 

Pl3 

=  S  25S  21 

=  S25 

+ 

s27, 

P  24 

=  ^26’ 

P  25 

=  *^24*^26 

=  V4 

+ 

S26’ 

Pa 5 

=  S22  S2A 

=  S22 

+ 

s24, 

P  21 

=  V0V2 

=  Vo 

+ 

s22, 

(36) 


ei  S20  S2i  S22  S23  S24  S25  S26  S27  S20  1 


(  52 1  522  )  (  ^23  ^24  )  (  ^25  526  )  ‘S27  ' 


(37) 


The  implementation  of  equations  (36)  and  (37)  in  Lab  VIEW  is  shown  in  Figure  35. 


Figure  35.  Lab  VIEW  Schematics  of  Channel  2  Position  Bit  and  Even  Residue  Flag. 
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3. 


Channel  3 


DeMorgan’s  Theorem  was  applied  to  equations  (25)  and  (28)  to  obtain 


and 
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e3  S30  Vl  S32  S33  S34  +  S15  S36  S37  S38 


=  S 


30 


( 53 1 532  ( 533  534  ( 535  536  ( 537  538  .  (39) 


The  implementation  of  equations  (38)  and  (39)  in  Lab  VIEW  is  shown  in  Figure  36. 
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s34  s35 


Figure  36.  Lab  VIEW  Schematics  of  Channel  3  Position  Bit  and  Even  Residue  Flag. 

4.  Sub-Channel  Flags 

Equation  (30)  was  used  to  construct  the  Sub-Channel  flags  in  Figures  37  and  38 
using  XOR  gates.  They  are  used  to  check  if  the  position  bits  of  each  RSNS  residue 
require  bit  reversal,  depending  on  which  sub-channel  it  is  from. 


Input  e ! 

% 

Output 

sc0 

Input  e3 

N 

XOR2 

Figure  37.  Sub-Channel  0  Flag. 
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Input  ej 

-1TV 

Output 

sc, 

Input  e2 

— ► — ))  N 

XOR2 

Figure  38.  Sub-Channel  2  Flag. 


C.  CONDITIONAL  BIT  REVERSAL 

Each  inverted  position  bit  was  mapped  from  two  possible  states,  depending  on 
whether  it  was  inverted  or  not.  For  example,  the  position  bit  p2 oa  was  mapped  from  p2o  if 
it  was  not  inverted  and  mapped  from  p21  if  it  was  inverted.  Thus,  the  Boolean  expression 


for  position  bit  p2 0a  can  be  expressed  as:  p20a  =  p21SC2  +  p20  SC2  =  ( p21SC2  j  j  p2()  SC2  j . 

The  conditional  bit  reversal  for  each  position  bit  can  be  accomplished  by  the  use 
of  a  2-to-l  multiplexer,  shown  in  Figure  39,  using  the  sub-channel  flag  as  a  control 
signal.  The  desired  output  is  either  Input  1  or  Input  2,  depending  on  whether  the  control 
signal  is  asserted  or  not  asserted,  respectively.  This  procedure  can  be  applied  to  all  the 
other  position  bits. 


Figure  39.  2-to-l  Multiplexer  (After  [4]). 


Given  Figure  39  as  a  basic  building  block,  the  multiplexer  architecture  is 
extended  to  form  the  Channel  2  and  Channel  3  bit  reversal  circuits,  shown  in  the  next  two 
sections.  The  only  modifications  are  that  NAND  gates  are  predominantly  used  to  achieve 
the  same  logic  operation  as  the  AND  and  OR  gates  in  Figure  39. 
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1. 


Channel  2  Reversal 


Recall  that  the  position  bits  of  Channel  1  are  never  reversed  from  Figure  21.  Thus, 
a  Channel  1  reversal  circuit  is  not  required  for  this  implementation.  For  Channel  2,  the 
position  bits  are  reversed  if  the  residue  is  from  Sub-Channel  2,  i.e.,  when  SC2  is  asserted, 
from  a  comparison  of  Figures  23  and  25.  The  Lab  VIEW  schematics  of  the  Channel  2 
conditional  bit  reversal  circuit  are  shown  in  Figure  40. 
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Figure  40.  Lab  VIEW  Schematics  of  Channel  2  Conditional  Bit  Reversal. 


2.  Channel  3  Reversal 

For  Channel  3,  the  position  bits  are  reversed  if  the  residue  is  from  Sub-Channel  1 

or  2,  i.e.,  when  SC0  is  asserted,  from  a  comparison  of  Figures  27  and  29.  The  Lab  VIEW 

schematics  of  the  Channel  3  conditional  bit  reversal  circuit  are  shown  in  Figure  41. 
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Figure  41.  Lab  VIEW  Schematics  of  Channel  3  Conditional  Bit  Reversal. 

D.  ALIGNMENT  LOGIC 

The  Lab  VIEW  schematics  of  the  alignment  logic  circuit  to  achieve  the  full  DR  of 
133  are  shown  in  Figure  42.  It  is  a  direct  mapping  of  equation  (32).  The  only 
modification  is  the  insertion  of  inverters  to  allow  the  proper  functioning  of  the  encoder 
circuit. 
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Figure  42. 
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Lab  VIEW  Schematics  of  Alignment  Logic  (Full  DR). 


Note  that  the  last  two  LPS  equations  ( LPS21  and  LPS22)  are  not  required  for  the 
truncated  DR  case,  and  have  to  be  removed  if  it  is  implemented. 
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E.  ENCODER 

Since  only  one  of  the  LPSX  equations  will  be  asserted  at  any  one  time  in  the 
Alignment  Logic  circuit  in  Figure  42,  an  encoder  can  be  used  to  convert  the  23  possible 
outputs  to  a  5 -bit  representation. 

To  achieve  this,  a  logic  table  is  derived  in  Table  2.  The  function  of  the  encoder  is 
to  convert  the  LPS  number  of  the  asserted  NAND  gate  into  a  five-bit  output.  For 
example,  if  n22  is  the  active  NAND  gate,  the  encoder  output  will  be  10110. 


Table  2.  Encoder  Logic  Table. 
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The  output  bits  of  the  encoder  are  thus  the  OR  combination  of  all  the  logic  Is  in 
their  respective  columns  in  Table  2.  For  example,  g4  will  be  asserted  if  any  of  the  nl6, 
nl7,  nl8,  nl9,  n20,  n21  and  n22  gate  is  active.  Again,  note  that  the  last  two  gates  are  not 
required  for  the  truncated  DR  case  and  have  to  be  removed  if  it  is  implemented. 

The  Lab  VIEW  implementation  of  the  encoder  for  the  truncated  DR  and  the  full 
DR  cases  are  highlighted  in  the  next  two  sections. 

1.  Truncated  Dynamic  Range 

The  logical  binary  expressions  for  the  21-to-5  encoder  for  the  truncated  DR  case, 
derived  from  Table  2,  are  given  by 


go  =  (n,  +n,  +n5)(n7  +  n9  +  nu)(nn  +  nl5  )(nl7  +  nl9 ) 

8l  =  (  n2  +  «3  +  n6  )  (  ni  +  «10  +  «1 1  )  (  nH  +  «15  )  ( ' «18  +  «19  ) 

g2  =  (n4  +  n5  +n6)(n7  +  n]2  +  nu  )(nu  +  n15  +  n20 )  (40) 

g3  =  (ns  +  n9  +  nw  )(nu+  nl2  +  nn )  (nl4  +  n]5 ) 
g4=(nu+nll+nl8)(n19+n20). 

The  implementation  of  equation  (40)  in  Lab  VIEW  as  an  encoder  circuit  to 
achieve  the  truncated  DR  of  128  is  shown  in  Figure  43. 
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Lab  VIEW  Schematics  of  Encoder  (Truncated  DR). 


2.  Full  Dynamic  Range 

The  logical  binary  expressions  for  the  23-to-5  encoder  for  the  full  DR  case, 
derived  from  Table  2,  are  given  by 


g0  =(«,  +n3+n5)( 

n7  +n9  +nn){ 

nx3  +  nl5  +n17 )( 

>!9+«2l) 

gl=(n2+n,+n6)( 

n7+nl0  +  nu) 

(nu  +  n]5  +  «l8) 

(n19  +  n22) 

g2  =(n4+n5+n6)( 

'n7+n]2+n]3) 

(nl4+nl5+n20 

)(n2X+n22) 

gj  =(ns+n9+nH>) 

(nu+nl2+nl3 

){nl4+nl5) 

g-4  =(«16  +«17  +«18 
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The  implementation  of  equation  (41)  in  Lab  VIEW  as  an  encoder  circuit  is  shown 
in  Figure  44. 


Figure  44.  Lab  VIEW  Schematics  of  Encoder  (Full  DR). 
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F. 


ADDER 


The  LabVIEW  schematics  of  the  adder  circuit  to  convert  the  encoder’s  5-bit 
representation  to  a  binary  representation  of  the  position  h  within  the  DR  are  shown  in 
Figure  45.  It  is  a  direct  implementation  of  Figure  32. 


Figure  45.  LabVIEW  Schematics  of  Adder. 

The  implementation  of  the  logic  blocks  of  the  RSNS-to-binary  converter  in  LabVIEW 
was  shown  in  this  chapter.  The  simulation  and  testing  of  this  converter  to  verify  that  the 
conversion  logic  is  functioning  properly  is  documented  in  the  next  chapter. 
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V.  SIMULATION  OF  RSNS-TO-BINARY  CONVERSION 


The  simulation  and  testing  of  the  RSNS-to-binary  converter  to  verify  that  the 
conversion  logic  is  functioning  properly  is  contained  in  this  chapter.  The  first  section  of 
this  chapter  is  a  description  of  the  creation  of  thermometer  codes  for  each  channel  in 
Lab  VIEW  to  simulate  as  the  inputs  to  the  conversion  algorithm.  The  second  section 
contains  the  key  results  obtained  for  three  simulation  models  (truncated  DR,  full  DR,  and 
full  DR  with  LPS  Priority  cases). 

A.  LAB  VIEW  THERMOMETER  CODE  GENERATORS 

The  Lab  VIEW  software  does  not  have  functions  to  generate  the  thermometer 
codes  directly.  However,  the  one-to-one  correspondence  between  the  RNS  and  RSNS 
vectors  can  be  exploited  to  convert  a  RNS  structure  into  a  RSNS  thermometer  code. 

Previously,  it  was  demonstrated  in  Figure  14  that  each  RSNS  thermometer  code 
can  be  mapped  to  a  unique  RNS  state  such  that  there  is  no  ambiguity  within  a  single 
channel  period.  Thus,  Lab  VIEW  functions  can  be  created  to  act  as  counters  to  cycle 
through  the  RNS  states  for  each  channel,  and  digital  logic  can  be  added  to  map  each  RNS 
state  to  the  required  RSNS  thermometer  code.  The  next  three  sections  contain  the  details 
of  the  procedure  to  generate  the  thermometer  codes  for  each  of  the  three  channels. 

1.  Channel  1 

The  logic  table  for  Channel  1  is  shown  in  Table  3.  The  RNS  states  are  represented 
by  the  bits  {/,  3 ,  2 ,  /, , ,  /,  0 }  and  are  the  inputs  to  the  Channel  1  thermometer  code  generator. 
The  desired  thermometer  code  output  is  represented  by  the  RSNS  bits 

{^16  ’  ^15  ’  ^14 ’^13 ’^12 ’‘hi’ ^10  }  ‘ 

Karnaugh  maps  are  used  to  carry  out  logic  minimization  for  Table  3,  shown  in 
Figure  46. 
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Table  3. 


Channel  1  (mx  =  7  )  Logic  Table. 
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Channel  1  Logic  Minimization  of  Table  3  (After  [4]). 


Figure  46. 
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The  logic  minimization  in  Figure  46  is  used  to  produce  the  logic  equations 


s 

s 

s 

s 

s 

s 

s 


10 

11 

12 

13 

14 

15 

16 


Z3  I  Z2  I  Zj  I  Zq  z3  I  z2  1  Zj  I  Zq 


h  44  4  4  4 4  zi  I  hl 2  1 1  4 4  1 1  6 4 


'  44  "*"44  "*"44  I  Z3Z2  )  1  Z3 Z2 


(Vo) 


44  "*~  4  4  4  "*"  4  4  4  I  44  11444  1 1  4  4  4 


44  "*"  444  "*"  4  4  4  (*26  )  ( 444  )  ( 4  4  4 


Z2  Zj  I  z3  Z2  Z- 


i4  (44) 


■  Z2Z^Zq  z2  H-  Zj  “l-  Zq 


(42) 


These  logic  equations  are  used  to  implement  the  Channel  1  thermometer  code 
generator  in  Lab  VIEW,  shown  in  Figure  47. 


Figure  47.  Lab  VIEW  Channel  1  Thermometer  Code  Generator. 


61 


Note  that  DeMorgan’s  Theorem  was  used  to  achieve  a  convenient  pipelined 
implementation  of  all  bits  in  the  three  channels. 


2.  Channel  2 

The  logic  table  for  Channel  2  is  shown  in  Table  4.  The  RNS  states  are  represented 
by  the  bits  \i2-,,i22,i2ri2()  }  and  are  the  inputs  to  the  Channel  2  thermometer  code  generator. 
The  desired  thermometer  code  output  is  represented  by  the  RSNS  bits 

{^27  ’  S26  ’  ^25  ’  ^24  ’  ^23  ’  ^22  ’  ^21  ’  ^20  }  • 


Table  4.  Channel  2  ( m2  =  8  )  Logic  Table. 
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Karnaugh  maps  are  used  to  carry  out  logic  minimization  for  Table  4,  shown  in 
Figure  48. 
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Figure  48.  Channel  2  Logic  Minimization  of  Table  4  (After  [4]). 
The  logic  minimization  in  Figure  48  is  used  to  produce  the  logic  equations 


S 2q  I  ^2  I  I  /q  I  I  I  /q 


52i  Z3Z2  +  +  h?o  I  *3*2  II  *3*1  11  *2*1  II  *i*o 


S22  h  h  L  *i  *3  Vo  *3*2*1  (*3*2  II  *2*1  I  *3*1*0  I  *3*2*1 


S23  l3l2+l3l2+l3l^  1  *3*2  II  *3*2  II  *3*1*0 


S24  l3 12  +  ?3?2?o  +  hhh  1  *3  *2  I  *3*2*0  I  *3*2*1 


S25  *3*2*1  *3*2*1  *3*2*0  (*3*2*1  (*3*2*1  (*3*2*0 


526  ?3  ?2  ?3^2^1?0  I  *3*2*1  II  *3*2  Vo 


(43) 


^27  i0  ^3  ^  ^  I 
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These  logic  equations  are  used  to  implement  the  Channel  2  thermometer  code 
generator  in  Lab  VIEW,  shown  in  Figure  49. 


Figure  49.  Lab  VIEW  Channel  2  Thennometer  Code  Generator. 
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3. 


Channel  3 


The  logic  table  for  Channel  3  is  shown  in  Table  5.  The  RNS  states  are  represented 
by  the  bits  {i34,i33,i32,in,i30}  and  are  the  inputs  to  the  Channel  3  thermometer  code 
generator.  The  desired  thermometer  code  output  is  represented  by  the  RSNS  bits 

{^38  ■>  *37  ■>  *36  ’  *35  ’  *34  ’  *33  ’  *32’  *31’  *30 J  • 

Table  5.  Channel  3  (m3  =  9)  Logic  Table. 
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Karnaugh  maps  are  used  to  carry  out  logic  minimization  for  Table  5,  shown  in 
Figure  50. 
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Channel  3  Logic  Minimization  of  Table  5  (After  [4]). 
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Figure  50. 


The  logic  minimization  in  Figure  50  is  used  to  produce  the  logic  equations 


530  =  Z4  +  Z3  + 12  + 1 


h  +io  ~  (h+h+h)(h  +io/ 


s3l=i4i0+i3+i2+i1  =  i4i0  *3+i2+*i 


) 


iS*32  ^3  I  ^  ^ 


do  (i3 +4j(Vo^ 


S33  i3i2+hh+hh+hh  I  hh  II  Mi  II  Z3Z2  II  z3zo 


s 34  hh+hh  +  hhh+hhh  I  hh  II  hh  II  z2zdo  II  hhh 


(44) 


^35  hh  +  hhlo+hhli  I  hh  II  z3zizo  II  hhh 


^36  hhhh+hh  I  z3z2zdo  II Z3Z2 


S37  hhh+hhh  I  hhl\  1 1  hhh> 


Sis  hhhh  h  h  “*~zi  “*~zo 


These  logic  equations  are  used  to  implement  the  Channel  3  thermometer  code 
generator  in  Lab  VIEW,  shown  in  Figure  5 1 . 
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Figure  51.  Lab  VIEW  Channel  3  Thennometer  Code  Generator. 
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B. 


SIMULATION  MODEL 


To  test  the  RSNS-to-binary  conversion,  a  simulation  model  was  created  with  the 
thennometer  code  generators  from  the  three  channels  enclosed  within  a  case  structure, 
shown  in  Figure  52.  A  ‘Master  Clock’  global  variable  was  created  to  determine  the 
starting  channel  and  to  step  through  each  of  the  three  channels.  Three  ‘Channel  Count’ 
global  variables  were  used  to  set  the  channel  start  positions. 


Ichannel  Start  Positions! 


flTf  |$Ch  1  Count 


[c^| —  |  $Chi  2  CounT 


[p~| —  |  ^>Chi  3  Count 


I  Startling  Channel  1 


fil — |  ^Master  Clock  | 


Figure  52.  Lab  VIEW  Simulation  of  RSNS-to-B  inary  Conversion. 
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c. 


SIMULATION  RESULTS 


After  running  the  simulation  model  in  Figure  52,  the  thermometer  code  inputs  and 
the  RSNS-to-binary  conversion  output  can  be  viewed  on  the  respective  scopes  located  in 
the  Lab  VIEW  Front  Panel. 

1.  Truncated  Dynamic  Range 

The  thennometer  code  inputs  and  the  simulated  RSNS-to-binary  output  with  no 
ambiguities  within  the  truncated  DR  of  126  are  shown  in  Figure  53.  The  RSNS  position 
index  h  is  from  733  to  857. 


Figure  53.  Simulated  RSNS-to-B inary  Output  (Truncated  DR  Case). 
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2. 


Full  Dynamic  Range 


To  achieve  the  full  DR,  the  encoder  circuit  in  Figure  44  was  used  in  place  of 
Figure  43  in  the  RSNS-to-binary  conversion.  The  ambiguities  in  the  simulated  RSNS-to- 
binary  output  within  the  full  DR  of  133  are  shown  in  Figure  54.  The  RSNS  position  index 
h  is  from  733  to  865. 


Figure  54.  Simulated  RSNS-to-Binary  Output  (Full  DR  Case  with  Ambiguities). 


3.  LPS  Priority  Circuit 

The  ambiguities  in  Figure  54  occur  because  the  DR  does  not  start  with  a  Sub- 
Channel  0  vector  and  end  with  a  Sub-Channel  2  vector,  as  required  in  Figure  9.  As  a 
result,  a  start  position  with  a  Sub-Channel  0  vector  (/?,  =  732)  has  to  be  chosen  outside 
the  DR,  during  the  DR  compression  process  in  Chapter  II. 
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As  the  DR  is  not  truncated  to  a  length  evenly  divisible  by  six,  there  may  be 
ambiguities  where  two  LPS  equations  will  be  asserted  at  the  same  time  -  one  of  the 
boundary  LPS  equations  (i.e.  LPSO  or  LPS22 )  and  LPSX.  If  the  full  DR  is  to  be 
maintained  with  no  ambiguities,  the  boundary  LPS  equation  should  always  be  ignored  in 
favor  of  the  LPSX  solution. 


This  can  be  done  by  implementing  an  encoder  with  an  addition  of  a  LPS  Priority 
circuit  to  check  for  such  cases,  and  giving  priority  to  the  center  2 1  NAND  gates  over  the 
first  and  last  NAND  gates.  This  is  shown  in  Figure  55. 


Figure  55. 


Lab  VIEW  Schematics  of  Encoder  with  LPS  Priority  Circuit. 
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4.  Full  Dynamic  Range  with  LPS  Priority  Circuit 

To  achieve  the  full  DR  with  no  ambiguities,  the  encoder  circuit  in  Figure  55  was 
used  in  place  of  Figures  43  and  44  in  the  RSNS-to-binary  conversion.  From  Figure  56,  it 
can  be  seen  that  there  are  no  ambiguities  in  the  simulated  RSNS-to-binary  output  within 
the  full  DR  of  133.  The  RSNS  position  index  h  is  from  733  to  865. 


Figure  56.  Simulated  RSNS-to-Binary  Output  (Full  DR  Case  with  No  Ambiguities). 

Table  6.  Comparison  of  Dynamic  Range. 


Model 

Position  Index  h 

Dynamic  Range 

Ambiguities? 

Truncated  DR 

733  -  857 

126 

No 

Full  DR 

733  -  865 

133 

Yes 

Full  DR  with  LPS  Priority 

733  -  865 

133 

No 
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A  comparison  of  the  dynamic  range  achieved  by  the  three  different  cases  is  shown 
in  Table  6.  The  simulation  proves  that  an  eight-bit  DR  of  133  can  be  achieved  for  a  three- 
channel  RSNS  ADC  with  Moduli  m\  =  7,  mi  =  8  and  m3  =  9.  This  is  in  agreement  with 
the  theoretical  DR  calculated  from  (5),  and  a  one -bit  improvement  over  that  achieved  in 
[4]- 

Simulation  of  the  RSNS-to-binary  conversion  algorithm  was  successfully  carried 
out  to  verify  that  it  is  working  properly  in  this  chapter.  In  the  next  chapter,  the 
implementation  of  the  comparator  circuit  and  RSNS-to-binary  conversion  algorithm  on  a 
FPGA  to  allow  the  code  to  run  at  a  higher  sampling  rate  and  frequency  is  elaborated.  The 
FPGA  is  then  integrated  with  the  front-end  photonics  implementation  to  form  the  overall 
folding  ADC  architecture. 
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VI.  ADC  INTEGRATION 


The  integration  process  to  form  an  overall  folding  ADC  system  is  documented  in 
this  chapter.  The  first  section  of  this  chapter  contains  an  illustration  of  the 
implementation  of  the  comparator  circuit  design  and  RSNS-to-binary  conversion 
algorithm  on  an  NI  FPGA  to  form  the  DDS  module.  The  integration  of  the  DDS  module 
with  the  front-end  PES  module  in  [5]  is  highlighted  in  the  second  section. 

A.  DIGITAL  DECODING  SUB-SYSTEM 

An  overview  of  the  DDS  module  from  the  bank  of  comparators  to  the  FPGA  is 
shown  in  Figure  57.  An  NI-9215  Analog  Input  Module  is  used  to  connect  the  three 
analog  photo-detector  inputs  from  the  front-end  PES  module  to  the  NI-9  111  Chassis  with 
an  onboard  XILINX  Virtex-5  LX30  FPGA.  The  FPGA  output  is  then  sent  to  the  NI-9012 
CompactRIO  (cRIO)  Real-Time  Controller  (RTC),  and  can  be  saved  to  a  file  using  a 
Lab  VIEW  Host  Interface  VI. 


Figure  57.  DDS  Module  Setup. 
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The  comparator  speed  was  a  limiting  factor  in  the  ADC  system  speed  as  actual 
comparator  ICs  were  used  for  sampling  outside  the  FPGA  in  [3].  The  key  improvement 
made  in  the  DDS  module  setup  in  Figure  57  is  that  the  comparator  circuits  reside  in  the 
FPGA,  allowing  them  to  sample  at  a  rate  equal  to  the  FPGA  speed. 

Lastly,  the  comparator  circuit  and  RSNS-to-binary  conversion  logics  are  run  on 
the  FPGA  to  guarantee  a  higher  FPGA  execution  speed  as  opposed  to  running  it  on  the 
RTC  with  a  lower  processing  speed. 

B.  COMPARATORS 

In  the  actual  ADC  architecture,  comparators  are  used  to  convert  the  analog  photo¬ 
detector  outputs  into  thermometer  codes  for  each  channel  in  place  of  thermometer  code 
generators.  The  threshold  values  for  the  comparators  can  be  expressed  as  [3]: 


T(k,mi)  =  Vfii  *  cos2 


K 

2 


/  7i  (kN  +  1)  n 

P  ?  P 

V  RSNS  RSNS  J 


(45) 


where  k  =  0,  1,  1,  VfS  =  full-scale  voltage  for  channel  Moduli  m„  N=  number  of 

channels,  and  Rrsns  =  period  of  RSNS  channel  =  2Nm<. 

The  full-scale  voltage  is  defined  as  the  maximum  amplitude  of  the  modulated 
signal  for  each  channel  and  detennined  in  [5]. 

1.  Channel  1 

Using  (45)  with  Vjs  =  2.5  (from  [5]),  we  calculate  the  quantized  threshold  values 
for  Channel  1  as: 


T{k,mx)  =  [ 0.0313,  0.2727,  0.7076,  1.2500,  1.7924,  2.2273,  2.4687],  (46) 
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The  implementation  of  the  Channel  1  comparators  in  LabVIEW  using  the 
threshold  voltages  calculated  in  (46)  is  shown  in  Figure  58. 


Figure  58.  LabVIEW  Schematics  of  Channel  1  Comparators. 


2.  Channel  2 

Using  (45)  with  V ^  =  5. 1 1  (provided  from  [5]),  we  calculate  the  quantized 
threshold  values  for  Channel  2  as 


T(k,m2)  =  [ 0.0419,  0.4306,  1.1355,2.0565,3.0535,3.9745,4.6794,5.0609],  (47) 


The  implementation  of  the  Channel  2  comparators  in  LabVIEW  using  the 
threshold  voltages  calculated  in  (47)  is  shown  in  Figure  59. 
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Figure  59.  Lab  VIEW  Schematics  of  Channel  2  Comparators. 

3.  Channel  3 

Using  (45)  with  =  4.96  (provided  from  [5]),  we  calculate  the  quantized 
threshold  values  for  Channel  3  as 

T(k,m 3)  =  [0.0377,  0.3323,  0.8859,  1.6318,2.4800,3.3282,4.0741,4.6277,4.9223],  (48) 

The  implementation  of  the  Channel  3  comparators  in  LabVIEW  using  the 
threshold  voltages  calculated  in  (48)  is  shown  in  Figure  60. 
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Figure  60.  Lab  VIEW  Schematics  of  Channel  3  Comparators. 

C.  FPGA  IMPLEMENTATION 

In  order  to  implement  the  comparators  and  RSNS-to-binary  conversion  on  the 
FPGA,  a  Lab  VIEW  FPGA  project  was  created  as  ‘FPGA789.1vproj’.  An  overview  of  this 
FPGA  project  is  shown  in  Figure  61. 

The  FPGA  has  a  Dynamic  Memory  Allocation  (DMA)  First-In-First-Out  (FIFO) 
buffer,  which  is  used  to  transfer  data  from  the  FPGA  to  a  host  computer  for  this  project. 
Lab  VIEW  also  requires  the  creation  of  an  FPGA  VI  to  run  the  necessary  codes  on  the 
FPGA  and  a  Host  VI  to  communicate  with  the  FPGA  VI,  which  will  be  explained  in  the 
next  two  sections. 
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Figure  6 1 .  Lab  VIEW  FPGA  Project  Overview. 


1.  FPGA  VI 

In  the  FPGA  VI,  the  three  photo-detector  inputs  were  read  in  using  the  NI-92 1 5 
Analog  Input  Module,  sampled  using  the  comparator  circuits,  and  processed  using  the 
RSNS-to-binary  conversion  algorithm.  The  RSNS-to-binary  output  was  then  stored  in  the 
DMA  FIFO.  The  number  of  elements  in  the  FIFO  can  be  set  in  Figure  61  by  double¬ 
clicking  on  the  FIFO  icon  and  changing  its  properties.  The  Lab  VIEW  schematics  of  the 
FPGA  VI  are  shown  in  Figure  62. 

A  point  to  note  is  that  the  thermometer  code  displays  were  removed  from  the 
Front  Panel  of  the  VI,  as  it  slowed  down  the  FPGA  execution  speed  when  it  was  required 
to  display  any  data. 
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The  truncated  DR  case  was  used  for  the  RSNS-to-binary  conversion  for  this 
integration  as  the  interferometers  used  in  [5]  could  not  create  enough  folds  of  the 
modulated  signals  to  exploit  the  whole  DR. 


Figure  62.  Lab  VIEW  Schematics  of  FPGA  VI. 
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2. 


Host  Interface  VI 


In  the  Host  VI,  Lab  VIEW  requires  that  a  reference  to  the  FPGA  VI  be  opened 
before  the  Host  VI  can  communicate  with  the  FPGA  VI.  Thereafter,  a  ‘FIFO  Read’ 
action  was  invoked  by  the  Host  VI  to  read  the  DMA  FIFO  data  in  the  FPGA  VI.  The  data 
was  then  saved  in  a  text  file,  using  the  ‘Write  to  Measurement  File’  VI.  The  Lab  VIEW 
schematics  of  the  Host  VI  are  shown  in  Figure  63. 


Figure  63.  Lab  VIEW  Schematics  of  Host  VI. 

3.  Results 

After  compiling  and  running  the  FPGA  and  Host  Vis,  proper  alignment  of  the 
modulated  signals  had  to  be  carried  out  to  ensure  that  the  RSNS  vectors  were  lined  up 
correctly  to  achieve  accurate  decoding.  The  DDS  outputs  for  a  1-kHz  triangular  and  1- 
kHz  sine  input  signal  are  shown  in  Figures  64  and  65,  respectively. 
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Time 

Figure  64.  DDS  Output  for  a  1-kHz  Triangular  Input  Signal. 


Figure  65.  DDS  Output  for  a  1-kHz  Sine  Input  Signal. 
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It  is  seen  from  Figures  64  and  65  that  the  DDS  module  was  able  to  decode  the 
input  signals  with  some  signal  clipping  and  uneven  step  sizes  in  the  output  waveforms. 

Implementations  of  the  comparator  circuit  and  RSNS-to-binary  conversion 
algorithm  on  the  FPGA  to  form  the  DDS  module  and  DDS  module  integration  with  the 
front-end  PES  module  to  form  the  overall  folding  ADC  architecture  were  carried  out  in 
this  chapter.  An  analysis  of  the  ADC  performance  characteristics  is  provided  in  the  next 
chapter. 
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VII.  ADC  PERFORMANCE 


In  this  chapter,  a  number  of  important  parameters  that  describe  an  ADC’s 
performance  are  presented.  Differential  and  integral  linearity  errors  are  plotted  for  the 
ADC  to  analyze  the  linearity  errors.  The  ADC  dynamic  range  is  determined  using  a  full- 
scale  sinusoid  and  a  Fourier  spectrum  analysis  of  the  noise  floor.  Finally,  the  ADC 
dynamic  performance  is  characterized  by  the  Signal-to-Noise  Ratio  (SNR),  SNR  plus 
distortion  (SINAD),  Total  Harmonic  Distortion  (THD)  and  the  Effective  Number  of  Bits 
(ENOB)  parameters. 

A.  LINEARITY  ERRORS 

A  comparison  of  the  ADC  1-kHz  sine  input  signal  and  the  DDS  output  signal  is 
shown  in  Figure  66,  where  it  can  be  seen  that  the  DDS  output  signal  is  able  to  follow  the 
input  signal,  with  some  quantization  and  linearity  errors. 


PES 1  kHz  Input  Signal  vs  DDS  Output  Signal 
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The  characteristic  transfer  function  of  a  1-kHz  triangular  waveform  is  used  to 
quantify  the  linearity  errors.  A  plot  of  the  ADC  transfer  function  using  a  1-kHz  triangular 
input  signal  is  shown  in  Figure  67.  The  quantization  errors  of  the  ADC  transfer  function 
in  Figure  67  are  shown  in  Figure  68.  Quantization  error  is  present  as  there  is  no  one-to- 
one  correspondence  between  the  input  and  output  voltage. 


Figure  68.  Quantization  Errors. 
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1. 


ADC  Resolution 


The  LSB  size  for  this  RSNS  ADC  is  defined  as: 

(49) 

M 

where  LSBrsns  is  the  ADC  resolution,  V/s  is  the  full-scale  voltage  and  M  is  the 
maximum  system  dynamic  range. 

2.  Differential  Non-Linearity 

The  Differential  Non-Linearity  (DNL)  is  expressed  as  [10]: 

DNLk  ={Vk-Vk_k)-LSB  (50) 

where  Vk  and  V/c_:  are  two  consecutive  code  transition  points,  and  ( Vk  -  Vk  , )  is  the  step- 

size.  The  DNL  is  the  maximum  deviation  in  the  output  step-size  from  the  ideal  value  of 
one  LSB. 


3.  Integral  Non-Linearity 

The  Integral  Non-Linearity  (INL)  is  defined  as  [10]: 

INLj  =  ^FJDNLk  =  Vj  -jLSB  (51) 

k= 1 

j 

where  V  is  the  sum  of  the  step-size  from  zero  to  the  jth  transition  point, 

k= 1 

and  jLSB  is  the  ideal  value  at  that  transition  point.  The  INL  is  the  maximum  deviation  of 
the  input/output  characteristic  from  a  straight  line  passed  through  its  end  points.  A  good 
ADC  typically  has  linearity  error  <  0.5  LSB.  [10] 

The  linearity  errors  (step-size,  DNL  and  INL)  of  the  quantized  signal  in  Figure  67 
are  shown  in  Figure  69.  The  step-size  plot  depicts  the  length  of  input  voltage 
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corresponding  to  each  quantization  level  as  the  input  voltage  increases.  It  shows  that  the 
maximum  step-size  is  2.7  LSB,  which  corresponds  to  the  code  transition  points  with  a 
maximum  DNL  value  of  1.7  LSB. 


A  maximum  INL  value  of  7.8  LSB  is  also  obtained  in  Figure  69.  The  INL 
increases  with  a  larger  input  voltage,  indicating  that  the  slope  of  the  ADC  transfer 
function  is  deviating  from  that  of  an  ideal  ADC  transfer  function  at  larger  input  voltages. 


Figure  69.  Linearity  Parameters  -  Step-Size,  DNL  and  INL. 
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B. 


NOISE  FLOOR  ANALYSIS 


The  process  for  computing  the  frequency  spectrum  and  analyzing  the  ADC  noise 
floor  is  illustrated  in  Figure  70. 

Magnitude 


Noise  Floor 
Harmonics 


Figure  70.  Process  for  Examining  ADC  Noise  Floor  (From  [10]). 

This  process  involves  passing  a  sinusoid  input  signal  through  the  ADC.  During 
quantization,  a  white  noise  process  (quantization  noise)  yk  is  added  to  the  analog  signal. 
The  ADC  noise  floor  is  analyzed  by  windowing  20  sets  of  4,096  digitized  samples  using 
a  Rlackman-Fl arris  window,  with  the  window  samples  represented  by  wk .  The 

Blackman-Harris  window  is  chosen  because  of  the  low  side-lobe  levels  it  can  achieve 

[10]. 

The  average  magnitude  spectrum  response  is  obtained  by  using  the  Discrete 
Fourier  Transform  (DFT)  to  transform  the  time-domain  signals  into  the  frequency 
domain  and  performing  asynchronous  point-by-point  spectral  averaging  of  the  20  sets  of 
data. 

The  purpose  of  conducting  this  noise  floor  analysis  is  to  determine  the  dominant 
noise  sources  in  the  ADC  system,  explained  in  the  next  two  sections. 

1.  Quantization  Noise 

The  theoretical  noise  floor  can  be  evaluated  by  examining  the  SNR  and 
considering  the  presence  of  quantization  noise  only.  By  normalizing  the  square  of  the 
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magnitude  of  the  spectral  average  with  the  fundamental  signal,  the  equation  to  calculate 
the  noise  floor  using  a  Blackman-Harris  window  is  [10]: 


F  2  =101og 


r3M^ 


10 


j 


+  6.02h  [dB] 


(52) 


where  n  =  Number  of  ADC  Bits  =  logi  M  =  logi  (41)  =  5.26  bits,  EB  =  Equivalent  Noise 
Bandwidth  of  Blackman-Harris  Window  =  2.0,  and  M=  Number  of  Samples  =  4096. 

The  theoretical  noise  floor  is  calculated  as  FM2  =  -64.13  dB  using  (52).  The  procedure  in 

Figure  70  is  then  used  to  obtain  actual  ADC  noise  floor  measurements  for  comparison. 

For  this  analysis,  a  1-kHz  sinusoidal  signal  was  sampled  at  100-kHz.  Twenty  sets 
of  4,096  digitized  signals  were  acquired  asynchronously  and  windowed  using  a 
Blackman-Harris  window.  The  MATLAB  Fast  Fourier  Transform  (FFT)  function  was 
used  to  compute  the  signal  spectrum.  The  point-by-point  spectral  average  of  the  twenty 
sets  of  data  was  then  calculated. 

The  spectral  average  of  the  1-kHz  sinusoidal  signal  is  shown  in  Figure  71,  using  a 
Blackman-Harris  window. 


Normalized  Magnitude  Squared  Spectrum  Using  Blackman-Harris  Window 

20, - ! - 1 - 1 - 1 - 1 - 1 - 1 - 1 - ! - 


Frequency(kHz) 

Spectral  Average  of  a  1-kHz  Sinusoidal  Signal. 
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Figure  71. 


Since  the  noise  floor  in  Figure  71  is  -64.5  dB,  which  is  close  to  the  theoretical 
noise  floor  of -64. 13  dB,  other  additive  noise  sources  (such  as  thermal  noise)  are  less 
dominant  in  the  system  compared  to  quantization  noise. 

2.  Clock  Jitter 

Two  frequencies  can  be  used  to  test  if  clock  jitter  is  a  dominant  noise  source,  with 
the  higher  frequency  being  twice  that  of  the  first  frequency  [10].  A  2-kHz  sinusoidal 
signal  was  sampled  to  compare  its  noise  floor  with  that  obtained  for  the  1-kHz  sinusoidal 
signal  in  Figure  71.  The  same  process  in  Figure  70  was  adhered  to  in  calculating  the 
magnitude  square  spectrum  of  both  signals. 

The  spectral  average  of  the  2-kHz  sinusoidal  signal  with  a  noise  floor  of  -59.5 
dB,  using  a  Blackman-Harris  window,  is  shown  in  Figure  72.  Since  the  difference 
between  both  noise  floor  levels  is  5  dB  (less  than  6  dB),  clock  jitter  is  not  the  dominant 
noise  source  [10].  This  is  expected  as  clock  jitter  is  not  expected  to  be  significant  at  these 
relatively  low  frequencies. 


Normalized  Magnitude  Squared  Spectrum  Using  Blackman-Harris  Window 


Spectral  Average  of  a  2-kHz  Sinusoidal  Signal. 
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Figure  72. 


c. 


DYNAMIC  PERFORMANCE  ANALYSIS 


Lastly,  the  ADC  performance  is  characterized  using  the  following  dynamic 
performance  parameters:  SNR,  THD,  SINAD  and  ENOB.  Two  different  frequencies  (1- 
kHz  and  2-kHz)  were  used  for  comparison. 

1.  Signal-to-Noise  Ratio 

The  ideal  SNR  equation  assumes  only  quantization  noise  and  is  expressed  as  [10]: 

SNR(dB)  -  6.02/7  +  1.76  (53) 

where  n  =  Number  of  ADC  Bits  =  log.,  M  =  log,  (41)  =  5.36  bits. 


2.  Total  Harmonic  Distortion 

The  THD  measures  the  hannonics  of  the  input  signal  that  show  up  at  integral 
multiples  of  the  fundamental  frequency,  and  is  a  measure  of  the  ADC’s  non-linearity.  It  is 
defined  as  [10]: 

THD(-dB)  =  20  log 

For  the  purpose  of  this  analysis,  the  input  signal’s  first  five  harmonics  were  measured  in 
decibels  and  used  in  (54). 


3.  SNR  Plus  Distortion 

The  SINAD  equation  takes  into  account  all  of  the  noise  (including  hannonics)  to 
give  an  indication  of  the  useful  ADC  dynamic  range,  but  excludes  the  DC  component.  It 
is  expressed  as  [10]: 


(55) 


4.  Effective  Number  of  Bits 

The  ENOB  is  a  measure  of  the  usable  ADC  dynamic  range,  which  is  reduced  due 
to  noise.  It  is  defined  as  [10]: 
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ENOB  = 


SIN  AD  - 1 .76  +  20  log 


/  full  scale  amplitude 
v  actual  input  amplitude  j 


6.02 


(56) 


5.  Summary  of  ADC  Dynamic  Performance  Parameters 

A  summary  of  the  calculated  parameter  values,  using  (53)  to  (56)  is  shown  in 
Table  7. 

Table  7.  ADC  Dynamic  Performance  Parameters. 


Frequency 

Ideal  SNR  (dB) 

SINAD  (dB) 

THD  (dB) 

ENOB  (Bits) 

1  kHz 

34.22 

34.21 

-61.68 

5.39 

2  kHz 

34.22 

34.21 

-61.01 

5.39 

From  Table  7,  it  is  observed  that  the  THD  for  the  2-kHz  signal  is  0.67  dB  higher 
than  that  for  the  1-kHz  signal.  The  SINAD  and  ENOB  values  are  similar  for  both 
frequencies,  indicating  that  there  was  insignificant  signal  distortion  when  the  frequency  is 
increased  from  1-kHz  to  2-kHz. 

An  enabler  to  achieving  this  is  the  design  of  comparator  circuits  in  the  FPGA, 
which  has  significantly  reduced  the  distortion  caused  by  sampling  errors,  as  compared  to 
using  actual  comparator  ICs. 

D.  SUMMARY 

In  this  chapter,  the  ADC  perfonnance  was  described  in  three  ways.  First,  the 
ADC  transfer  function  was  examined  for  linearity  errors.  Second,  the  noise  floor  was 
analyzed  to  detennine  the  main  sources  of  noise  in  the  system.  Lastly,  the  dynamic 
performance  parameters  of  the  ADC  system  were  characterized.  Key  conclusions  and 
recommendations  for  future  research  will  be  provided  in  the  next  chapter. 
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VIII.  CONCLUSION 


The  key  conclusions  obtained  from  this  project  and  recommendations  for  future 
research  are  provided  in  this  chapter. 

A.  KEY  CONCLUSIONS 

The  goal  of  this  thesis  was  to  conduct  hardware  and  software  implementation  of 
the  DDS  module  of  the  folding  ADC  architecture  (for  Moduli  7,  8,  9)  from  the  bank  of 
comparators  to  the  RSNS-to-binary  conversion  within  the  FPGA  as  well  as  integration 
with  the  front-end  PES  module  of  this  ADC  design.  This  was  accomplished  via  several 
milestones  described  below. 

Firstly,  the  RSNS  dynamic  range  computation  algorithm  in  [6]  and  [7]  was 
verified  to  be  correct,  proving  that  an  eight-bit  DR  of  133  can  be  achieved  theoretically 
for  a  three-channel  RSNS  ADC  with  Moduli  m.\  =  7,  im  =  8  and  m3  =  9. 

Secondly,  the  RSNS-to-binary  algorithm  was  implemented  in  LabVIEW  and 
shown  to  achieve  the  DR  value  of  133,  which  is  in  agreement  with  [6]  and  [7],  and  a  one- 
bit  improvement  over  that  achieved  in  [4].  Design  of  thermometer  code  generator  circuits 
and  simulation  of  this  algorithm  were  carried  out  to  verify  that  it  is  working  properly 
before  connecting  to  actual  signals. 

Thirdly,  the  comparator  circuits  and  RSNS-to-binary  conversion  algorithm  were 
implemented  on  the  FPGA,  allowing  the  ADC  to  achieve  a  higher  sampling  frequency. 

Lastly,  the  DDS  module  was  integrated  with  the  front-end  PES  module  in  [5]  to 
form  a  folding  ADC  system  and  characterization  of  the  ADC  performance  was  carried 
out.  Analysis  of  the  results  attributed  the  dominant  noise  source  in  the  ADC  system  to 
quantization  noise,  with  the  ADC  remaining  resilient  to  errors  caused  by  other  additive 
noise  sources  and  comparator  sampling. 
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This  electro-optic  RSNS  ADC  system  has  been  demonstrated  to  work  and 
produces  a  seven-bit  output  with  relatively  simple  hardware  and  software.  Due  to  the 
reduced  number  of  hardware  and  software  components  and  the  large  bandwidth  of  the 
photonics  components,  the  energy  and  size  savings,  as  well  as  increase  in  speed,  make 
this  folding  ADC  design  appealing  for  defense  applications,  such  as  unmanned  systems, 
direction-finding  antenna  architectures  and  electronic  warfare  system-on-a-chip 
applications. 

B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

While  this  thesis  is  concerned  with  optimizing  the  DDS  module  implementation 
of  the  folding  ADC  system,  various  system  trade-offs  had  to  be  made  for  integration 
purposes.  Further  upgrades  and  research  can  be  carried  out  in  various  aspects  of  the 
system,  detailed  below. 

1.  Bandwidth  Upgrade 

The  bandwidth  of  this  ADC  system  is  currently  limited  by  the  NI-9215  Analog 
Input  Module’s  data  rate  of  100  kS/s.  The  highest  frequency  it  can  sample  is  50  kHz, 
based  on  the  Nyquist  criteria.  There  are  several  wider-bandwidth  NI  modules  available  to 
replace  this,  such  as  the  NI-5761  Digitizer  Adapter  Module,  which  has  a  data  rate  of  250 
MS/s  and  can  sample  frequencies  up  to  125  MHz  [11],  The  only  drawback  is  that  choices 
of  wide-bandwidth  modules  are  limited  due  to  lack  of  commercial  development  and 
require  impedance-matching  for  maximum  power  transfer. 

2.  FPGA  Upgrade 

The  speed  of  this  ADC  system  is  currently  limited  by  the  XILINX  Virtex-5  LX- 
30  FPGA  on  the  NI-9111  Chassis.  There  are  several  higher-capacity  NI  FPGA  modules 
available  to  replace  this,  such  as  the  FlexRIO  PXIe-7965R  module,  which  has  a  XILINX 
Virtex-5  SX-95T  FPGA  with  a  clock  speed  of  550  MHz  and  the  ability  to  handle  single- 
ended  Input/Output  (I/O)  up  to  800  Mbps  [12].  This  module  did  not  arrive  in  time  for  the 
project  due  to  procurement  delay,  but  can  be  easily  substituted  when  it  is  available. 
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There  are  other  high-speed  FPGAs  (up  to  1  GHz)  that  could  not  be  used  as  they 
do  not  interface  with  NI  Lab  VIEW.  Nevertheless,  it  is  envisaged  that  an  Application- 
Specific  Integrated  Circuit  (ASIC)  can  be  developed  for  testing  once  the  FPGA  circuit 
design  has  been  fixed.  This  will  remove  the  constraint  of  having  to  rely  only  on  the  NI 
programming  environment. 

The  envisaged  DDS  setup,  after  incorporating  the  component  upgrades,  is 


illustrated  in  Figure  73. 


FPGA  Output 


Figure  73.  Upgrades  to  DDS  Setup. 


3.  Higher-Moduli  Configurations 

The  modular  structure  of  the  RSNS-to-binary  converter  allows  the  ADC  system  to 
be  easily  scalable  to  higher-moduli  configurations  or  configurations  with  more  than  three 
channels.  This  will  increase  the  ADC  resolution  and  enable  the  ADC  system  to  achieve  a 
dynamic  range  that  is  better  than  the  current  eight-bit  DR  of  133. 
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APPENDIX.  GENERALIZED  CHINESE  REMAINDER 
THEOREM  PROCEDURE  TO  SOLVE  FOR  COA  SHIFTS 


The  COA  shift  hs\  is  the  least  positive  solution  to  the  following  sets  of 
congruence  equations: 

~~~  =  0  (mod  7)  - >  hsl  =  1  (mod  2 1) 

=  0  (mod  8)  - >  hA  =  1  (mod  24)  (57) 

^sl  =  0  (mod  9)  - »  hsl  =  -2  (mod  27) 


Equation  (57)  cannot  be  solved  directly  using  the  standardized  Chinese 
Remainder  Theorem  (CRT)  as  the  three  moduli  are  not  Pair-Wise  Relatively  Prime 
(PRP).  The  generalized  CRT  procedure  must  be  used.  The  first  step  is  to  break  each 
equation  into  its  constituent  equations  [10]: 


1  (mod  21) 


jhs]  =  1  (mod  7) 
[hsl  =  1  (mod  3) 


h 


s  1 


1  (mod  24) 


|  hsl  =  1  (mod  8) 
lAi  =  1  (mod  3) 


(58) 


[/f.i  =  “2  (mod 


hsl  =  -2  (mod  27)  < 


hsl  =  -2  (mod 


3) 

3) 


[hs j  =  -2  (mod  3) 


The  second  step  is  to  group  all  constituent  equations  with  the  same  moduli 
together,  and  solve  for  the  remaining  congruence  equations  using  the  standardized  CRT 
method: 

hsl  =  1  (mod  7) 

hsl  =  1  (mod  8)  (59) 

hsX  =  -2  (mod  3) 
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The  standardized  CRT  solution  is: 


K  =  Z 


M , 

—bflt 

m. 


where  M  =  n*-(7)(8)(3)  =  168 

mi  =  Moduli  /  =  [mt  m2  ra3]  =  [7  8  3] 
a;.  =  Residue  of  Moduli  i  =  a2  a3]  =  [l 


(60) 


The  coefficients  b\  can  be  found  via  the  following  procedure: 


Find  h  : 


(m  1 

( 168 

— ,ml 

- ’7 

l  7  J 

=  (24,7) 


24(l)  +  7(0)  =  24 
24(0)  +  7(1)  =  7 


} 


24 


24 


24 


W-(o) 

(O)-(l) 


24 
.  7  . 

7  ' 
3 


+  7 


(o)-(i) 


fl¬ 


ay 


+  7 


24 
7  . 

7 

3 


jy 


=  24(l)  +  7(-3)  =  3 
=  24(-2)  +  7(7)  =  l  !• 


Find  bn : 


M 

—  ,m2 

Vm2  y 


168 


=  (21,8) 


v  o  y 


2 1  (l)  +  8  (0)  =  21 

2 1  (0)  +  8  (l)  =  8 


! 


21 


21 


21 


21 


21 


(i)-(o) 

v 

(o)-(i)' 

(l)-(-l) 

(-l)-(2) 


21 

\ 

_~8~ 

-J 

8 

y 

_5_ 

y_ 

'  5 

\ 

vb 

-J 

+  8 


(o)-(i) 


+  8 


W-(-2) 


21 

¥ 

8 


V,l5  jy 


+  8 


+  8 


(—2)  —  (3) 


Vl3  jy 


(3)-(-5) 


=  2 1  (l)  +  8  (—2)  =  5 
=  21(-l)  +  8(3)  =  3 
=  2l(2)  +  8(-5)  =  2 
=  21  (—3)  +  8(8)  =  1  \  :.b 


=  3 

2 

-2 


=  2 

1 

1 

1 

-  3 
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P.  ,,  (m  )  fl68  ^  . 

Find  o3 :  — ,m3  =  - ,3  =(56,3) 

ym3  )  V  3  J 

56(l)  +  3(0)  =  56 

56(0)  +  3(l)  =  3  }  y  =18 


(l)-(0)' 

V 

56  w  r  ( 

f  +3  (0)-(l) 

.L  JyJ  L  v 

56  \ 

-3  JJ 

=  56(l)  +  3(-18)  =  2  | 

(0)-(l)' 

V 

4  )j  +  3L-(-18) 

'  3  " 

J-2_, 

=  56  (—1)  +  3(19)  =  1  j 

,  M,  M,  M  ],  x 
/isl  =  — Ojflj  H b2ct2  H o3a3  (modM) 

J 

hsl  =  [24(-2)(l)  +  2l(-3)(l)  +  56(-l)(-2)](mod  168)  =  l(mod  168) 


Similarly,  the  COA  shift  hs2  is  the  least  positive  solution  to  the  following  sets  of 
congruence  equations: 

=  0  (mod  7)  - »  hs2  =  2  (mod  2 1) 

=  0  (mod  8)  - >  hs2=- 1  (mod  24)  (61) 

~~~~  =  0  (mod  9)  - >  hs2  =  -1  (mod  27) 


Equation  (61)  cannot  be  solved  directly  using  the  standardized  CRT  as  the  three 
moduli  are  not  PRP.  The  generalized  CRT  procedure  must  be  used.  Again,  the  first  step  is 
to  break  each  equation  into  its  constituent  equations  [10]: 
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Ki=  2  (mod  2 1)  < 


K 2  =  2  (mod  7) 
hs2  =  2  (mod  3) 


hs2  =  -1  (mod  24) 


\hs2  =  -1  (mod  8) 
A 2  =  -1  (mod  3) 


hs2  =  -1  (mod  27) 


Ki  =  -1  (mod  3) 
hs 2  =  -1  (mod  3) 
hs2  =  -1  (mod  3) 


(62) 


The  second  step  is  to  group  all  constituent  equations  with  the  same  moduli 
together  and  solve  for  the  remaining  congruence  equations  using  the  standardized  CRT 
method: 


hs2  =  2  (mod  7) 
/zs2  =  -1  (mod  8) 
hs2  =  -1  (mod  3) 

The  standardized  CRT  solution  is: 


Ki  =  — 

i= i  mt 

where  M  =  n»,=(7)(8)(3)  =  168 

ms  =  Moduli  /'  =  [/«,  /n2  /«,  ]  =  [7 

at  =  Residue  of  Moduli  i  =  [<;/,  a2 


8  3] 

«3]  =  [2  “I  -1] 


(63) 


(64) 


The  coefficients  b{  are  the  same  as  for  hsl ,  as  the  moduli  configuration  is  the  same 
for  both  COA  shifts: 


A  = 


f  M  ,  M  ,  M  ,  A 
—  6,a,  H - b2a2  H - h3a3 


v/«i 


(modM) 


7 


.'A2  =[24(-2)(2)  +  21(-3)(-l)  +  56(-l)(-l)](mod  168)  =  23(mod  168) 


(65) 
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